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Objective 

The  objectives  of  this  research  were  threefold:  (1)  to  provide  the  methodology  needed  to 
quantify  the  two-dimensional  variations  in  bottom  topography  and  three-dimensional 
variations  in  subbottom  properties  observed  within  the  ARSRP  Central  Atlantic  study  area; 
(2)  to  apply  this  methodology  to  geological  data  sets  collected  in  the  study  area;  and  (3)  to 
use  the  resulting  geological  models  to  understand  the  bottom  reverberations  observed  at  long 
ranges  across  the  study  area  and  other  pans  of  the  deep  ocean. 

Background 

Previous  research  by  marine  morphologists,  including  the  work  done  by  the  P.I.  and  his 
colleagues  [Goff  and  Jordan,  1988,  1989a,b,  1990;  Goff,  Jordan,  Edwards,  and  Fomari, 
1991],  has  demonstrated  that  realistic  models  of  the  seafloor  can  be  formulated  by  combining 
deterministic  descriptions  of  large-scale  features  with  stochastic  parameterizadons  of  small- 
scale  features.  Most  of  the  analyses  done  to  date,  however,  have  concentrated  on  stochastic 
descriptions  of  the  hard-rock  topography  generated  by  seafloor  spreading  at  ridge  crests.  In 
particular,  little  attention  has  been  given  to  the  variations  of  topography  and  subbottom 
properties  caused  by  sedimentation  and  post-depositional  sediment  transport,  which  play  a 
major  role  in  controlling  the  subsurface  properties  of  the  seafloor,  including  elastic  and 
anelastic  properties,  porosity,  and  permeability.  A  major  focus  of  our  recent  research  has 
been  to  develop  methods  for  describing  sediment  distributions  on  a  rough  seafloor  and  to 
incorporate  these  descriptions  into  environmental  models. 

Approach 

We  first  prescribe  sediment  deposition  and  post-depositional  transport  by  a  simply 
parameterized  model  and  solve  the  forward  problem:  given  the  values  of  the  parameters, 
which  include  a  description  of  the  underlying  basement  topography,  predict  the  sediment 
distribution  and  observable  aspects  of  seafloor  structure.  The  observables  are  (1)  the 
bathymetry  of  the  sediment-water  interface  obtained  by  a  multibeam  swath-mapping  system, 
and  (2)  the  apparent  sediment  thickness  measured  from  a  single-channel  3.5-kHz 
echosounder  and/or  seismic  profiler.  We  then  set  up  and  solve  the  inverse  problem:  given 
these  observations  and  a  description  of  their  uncertainties,  estimate  the  m<^el  parameters. 
Finally,  we  evaluate  and  interpret  the  results;  if  the  fit  to  the  data  is  satisfactory,  we  proceed 
with  a  geological  interpretation  of  the  model;  if  not,  we  seek  modifications  to  our  pa- 
rameterizations  that  deliver  better  descriptions  of  the  data,  and  we  return  to  the  inversion 
step. 


EssuUs 

The  results  during  this  period  of  performance  period  are  described  in  the  following 
documents;  which  are  included  as  sections  1, 2  and  3  of  this  report: 

Section  1:  Webb,  Helen  F.  and  Thomas  H.  Jordan,  Quantifying  the  Distribution  and 
T  msport  of  Pelagic  Sediments  on  Young  Abyssal  Hills,  submitted  to  Geophys. 
Res.  L£tt.  May  14, 1993. 


Section  2:  Webb,  Helen,  Modeling  Sedimentation  of  a  Rough  Seafloor,  submitted  as  partial 
fulfillment  of  the  requirements  of  the  Generals  Examination,  MIT/WHOI  Joint 
Program  in  Oceanography,  advised  by  Thomas  H.  Jordan,  August  5  1992. 

Section  3:  Goff,  John  A.,  Helen  F.  Webb  and  Thomas  H.  Jordan.  Stochastic  Interpolation  of 
Wide-Beam  Bathymetric  Profile  Within  the  ONR-ARSRP  Acoustic 
Reverberation  Corridor,  Office  of  Naval  Research  S.R.P.  on  Acoustic 
Reverberation,  Technical  Report  No.  2, 1992. 


Our  basic  forward-modeling  tools  are  efficient  numerical  algorithms  for  characterizing 
the  sediment  distribution  on  very  rough  (fractal)  models  of  seafloor  topography.  Sediment- 
particle  dynamics  are  described  by  a  diffusive  redistribution  process  in  which  the  diffusivity 
parameter  jf  and  the  sediment  flux  rate  F  can  vary  in  time  and  space.  We  have  developed  a 
3D  version  of  the  algorithm  and  implemented  it  on  a  massively  parallel  supercomputer,  the 
NRL  SOOO-processor  Connection  machine  [Webb  and  Jordan,  1992a].  Numerical 
simulations  have  been  run  for  various  values  of  k,  F,  and  the  five  parameters  of  a  Goff- 
Jordan  description  of  the  underlying  basement  topography.  From  a  stochastic  analysis  of 
these  simulations,  we  have  constructed  the  functional  relationships  between  the  model 
parameters  and  two  sets  of  observables:  (1)  the  apparent  RMS  height  H,  apparent 
characteristic  wavenumber  k ,  and  apparent  fractal  dimension  D  estimated  from  a  multibeam 
swath  by  the  inversion  algorithm  of  Goff  and  Jordan  [1988],  and  (2)  the  cumulative  slope 
distribution  function  P{s,Ax)  estimated  from  a  single-channel,  narrow-beam  echosounder 
sampled  with  a  horizontal  spacing  Ax . 

Using  these  functional  relationships,  we  have  inverted  the  data  collected  by  the  R/V 
Ewing  Hydrosweep  multibeam  system  during  the  ARSRP  Geology  and  Geophysics 
Reconnaissance  Cruise  in  Jul-Aug,  1992.  The  preliminary  map  of  the  average  sediment  load 
L  (the  time  integral  of  F)  for  the  ARSRP  survey  region  is  shown  in  Fig.  1 .  We  note  that  L  is 
highly  variable  across  the  study  area,  reflecting  spatial  variations  in  sea-floor  age,  variations 
in  F  controlled  by  temporal  variations  in  the  CCD,  and  lateral  transport  of  sediments 
controlled  by  large-scale  variations  in  bottom  topography. 

We  also  developed  a  DAT-based  system  for  the  digitization  and  playback  of  3.5-kHz 
profiler  returns.  This  system  was  installed  on  RtV  Ewing  and  successfully  collected  data 
continuously  during  the  ARSRP  Geology  and  Geophysics  Reconnaissance  Cruise.  These 
data  will  be  analyzed  for  small-scale  variations  in  reflectivity  associated  with  volumetric 
heterogeneities  in  sediment  ponds,  to  be  used  as  input  to  numerical  simulations  of  sediment 
backscattering  by  the  ARSRP  Modeling  Group. 
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FIGURE  1.  A  preliminary  map  of  the  average  sedimmit  load  L,  shown  by  the  numbers  in  meters,  obtained  by 
inverting  multibeam  slq)e-distribution  data  collected  for  the  ARSRP  survey  region  in  Jul-Aug,  1992.  L  is  highly 
varial^  across  the  study  area,  reflecting  spatial  variations  in  sea-floor  age,  variations  in  F  controlled  by  temporal 
variations  in  the  CCD,  and  lateral  tran^rt  of  sediments  crmcrolled  by  large-scale  variations  in  bottom  topography. 
The  data  fiom  some  areas,  including  the  large  "Site  A"  sediment  pond  (46.3**W,  26.3^,  have  not  been  analyzed. 


QUANTIFYING  THE  DISTRIBUTION  AND  TRANSPORT  OF  PELAGIC  SEDIMENTS 

ON  YOUNG  ABYSSAL  HILLS 

HELEN  F.  Webb  and  Thomas  H.  Jordan 
Dq>aitinent  of  Earth,  Atmospheric  and  Planetary  Sciences,  MIT,  Cambridge. 


Abstract.  In  pelagic  environments,  the  redistribution  of  sediments  can  be  modeled  as  a 
diffusive  process  in  which  material  is  transported  from  topographic  highs  to  lows  at  a  rate 
determined  by  an  apparent  diffusivity  k.  We  have  developed  a  statistical  methodology  for 
estimating  the  k  together  with  the  average  sediment  load  L  and  the  root-mean-square  height  H  qf 
the  abyssal-hill  basement.  Our  approach  is  to  generate  synthetic  basement  topognqphies  based  on 
the  Goff-Jordan  description  of  abyssal-hill  morphology,  distribute  sediments  over  these 
topographies  according  to  the  diffusion  assumption,  and  calculate  the  slope  distributions  of  the 
resulting  sediment-water  interfaces;  we  then  determine  x;  L,  and  //  by  fitting  the  slope-distribution 
function  to  the  values  observed  fnxn  narrow-beam  bathymetric  premies.  The  local  values  of  these 
parameters  have  been  derived  from  center-beam  Hydrosweep  data  collected  on  young  (S-25  Ma) 
abyssal  hills  in  the  Acoustic  Reverberation  Special  Research  Project  corridOT  on  the  west  flank  of 
the  Mid-Adantic  Ridge.  The  average  value  of  the  apparent  diffiisivi^  is  k  =  0.131.014  mVyr, 
which  is  the  first  determination  of  this  parameter  in  a  deep-sea  environment  H  varies  betweoi  188 
and  285  m,  with  a  mean  of  230  1  7  m.  L  is  estimated  to  be  as  thin  as  10  m  on  S-Ma  crast  near  the 
ridge  crest  thickening  to  values  of  greater  than  50  m  at  10-15  Ma  and  then  decreasing  to  less  than 
30  m  on  25-Ma  crust  in  the  western  portion  of  the  corridor.  The  decrease  in  sediment  load  with 
age  may  be  associated  with  a  shallower  carbonate  compensation  depth  prior  to  the  late  Miocene, 
hypothesized  by  previous  authors.  The  trends  in  our  sediment-load  estimates  agree  with  single¬ 
channel  data  collected  in  the  region. 
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Introduction 


Peb^  sedimentation  is  die  primary  modifier  of  seafloor  tqiogn^hy  on  the  flanks  of  the  mid¬ 
ocean  lidfes,  smoothing  and  eventually  burying  the  abyssal  hills  generated  by  seafloor  ^treading. 
Undentanding  the  distribution  of  sediments  <mi  young  seafloor  and  the  processes  that  govern  thdr 
po«-depositional  transport  are  central  problems  in  marine  morphology.  Thin  pelagic  sediments  are 
notoriously  difficult  to  study,  however,  because  the  scattering  by  rough  abyssal-hill  topography 
obscures  die  wide-beam  sediment  profiles  obtained  from  surface  ships.  In  this  report,  we  outline  a 
new  methoddogy  for  quantifying  sediment  distribution  and  transport  that  utilizes  cmly  bathymetric 
data  routinely  collected  by  narrow-beam  echosounders. 

Considered  in  detail,  die  processes  that  govern  deqy  sea  sedimentation,  such  as  the  generatitMi 
of  small-scale  turbidity  currents,  the  interactions  of  bottom  currents  with  pre-existing  topography, 
and  tnoturbation,  are  complex.  But  averaged  over  geological  time,  these  processes  act  to  transport 
sediments  from  topographic  highs  to  topographic  lows  at  all  scales  [e.g.  Marks.  1981].  We 
therefore  model  sedimentation  as  a  simple  diffusive  process  in  which  the  lateral  flux  of  sediments 
is  proportional  to  the  local  topographic  gradient  [Culling,  I960].  Diffusion  models  have  been  used 
to  eiqilain  the  evolution  of  scarp-like  landfOTms  [Hanks  et  al.,  1984],  as  well  as  the  stratigraphy 
found  in  frneland  basins  [Flemings  and  Jordan,  1989]  and  shallow  marine  settings  [Kaufman  et 
al.,  1991].  To  our  knowledge,  however,  sedimentatimi  in  pelagic  environments  has  not  been 
addressed  from  this  perspective. 


forward  Modeuno 

Let  h{x,t)  be  the  height  of  a  durable  basement  b{x)  covered  by  an  erodible  laver  of  sediment 
5(x,r),  let  F(x,r)  be  the  sedimentation  rate,  and  let  rr  be  the  diffusivity  of  lateral  sediment 
tranqwrt  The  sediment  transptnt  equatitm  is 

=  KVlhixj)+  F(x,t)  (1) 

where  is  the  surface  Lapacian.  Eqn.  (1)  is  to  be  solved  for  h(x,t)  =  d(x)  -i-  s(x,r)  subject  to 
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the  ooQStnint  tfiat  j(x,/)  ^  0.  The  inequality  makes  the  problem  highly  non-linear.  Ccnnplexity  is 
also  mtioduced  by  the  structure  of  the  basement  u^graphy  h(x) ,  which  has  a  fractal  character. 

We  have  obtained  representations  of  the  sediment-water  interface  by  numerically  integrating 
(1),  starting  from  basement  tqx)graphies  generated  as  realizations  of  the  stochastic  process  of  Gctff 
and  Jondan  [1988, 1990].  The  Gcrff-Joidan  nxrdel  represents  the  seafloor  as  a  spatially  stationary, 
anisoiR^ic  Gaussian  random  field  specified  by  five  parameters:  the  root-mean-square  (nns)  height 
die  characteristic  length  A,  and  width  of  the  abyssal  hills,  the  abyssal-hill  azimuth  and  a 
fractal  dimension  D.  For  the  applications  discussed  here,  we  fixed  the  scale  lengths  and  fractal 
dimension  at  values  typical  oi  the  western  flank  of  the  Mid-Adantic  Ridge  [Goff  and  Jordan,  1990; 
J.  Goff,  personal  communication,  1992]:  22  km,  A^  s  5.4  km,  and  D  =  2.2. 

The  computadons  were  performed  on  the  Naval  Research  Laboratory's  CM-2  Connection 
Madune  using  a  horizontal  grid  with  dimensions  of  100  km  x  100  km  and  a  grid  point  spacing  of 
Ax  =  100  m.  The  flux  rate  F  was  assumed  to  be  constant  in  space  and  time,  so  that  the  total 
sediment  load  for  a  seafloor  of  age  7  is  L  =  FT.  In  this  case,  the  diffusivity  scales  as  x  =  KLX/T, 
where  JIT  is  a  dimensionless  diffusivity.  Synthetic  topographies  wore  calculated  for  L  100  m  for 
a  wide  range  of  AT,  7,  and  H. 


Parameter  determination 

Two  examples  of  our  calculations  are  shown  in  Figure  1.  Low  values  of  ir  produce  seafloor 
morphologies  in  which  sediments  drape  hillsides,  while  higher  values  transport  the  sediments 
downslope,  exposing  the  peaks  and  creating  flat  sediment  ponds.  This  behavior  is  reflected  in  the 
slope  distribution,  which  is  a  sensitive  indicator  of  deviations  from  the  Gaussian  character  assumed 
fOT  the  basement  topography  [Shaw  and  Smith,  1990].  For  a  profile  of  constant  orientation  with 
respect  to  the  topographic  axis  Cs  sampled  at  a  constant  horizontal  interval  u,  we  define  the 
cumulative  slope  distribution  function  S{0,u)  as  the  fraction  of  profile  length  with  an  apparent 
slope  angle,  arctan[(h(x+u)-/i(x))/u],  less  than  0. 

Si0,u)  can  be  estimated  directly  from  single-channel,  narrow-beam  echograms  and,  as 
illustrated  in  Figure  2,  its  shape  is  diagnostic  of  //,  L,  and  k.  For  large  values  of  the 
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diineiisionless  diffusivity  and  for  sediment  loads  that  are  a  significant  fraction  of  the  basonent 
topography,  the  ponds  are  well  defined,  and  S(d.ii)  displays  a  sharp  Idnk  at  a  slqx  sngle  9^  that 
separates  the  pcmds  from  hillsides;  i.e..  5(0o.ti)  measures  die  fraction  of  the  profile  occupied 
ponded  sediments.  This  fraction  decreases  as  L  decreases  and  the  rate  of  this  decrease  is  most 
rqiid  when  UH  is  small.  The  slope  of  S{9m)  at  9  «  6*0  is  sensitive  to  x;  and  the  sl<^  of  S{9ju) 
at  &»  6*0  depends  primarily  tm  //.  In  other  wOTds.  by  least-squares  fitting  the  theOTetical  values 
of  5(9.m)  to  the  observed  values,  we  can  determine  //.  L.  and  k.  Although  the  results  are 
ctmditional  cm  the  values  assumed  for  A^.  and  f>.  which  we  have  fixed  in  our  inversions, 
numerical  experiments  show  that  5(d.u)  is  only  vwakly  dependent  on  these  additional  parameters. 

APPUCATIONTOTHE  ARSRP  CORRIDOR 

This  procedure  has  been  tested  on  a  data  set  collected  during  Cruise  9208  of  the  R/V  Maurice 
Ewing  (14  Jul-18  Aug  92)  as  part  of  the  Office  of  Naval  Research  (ONR)  Acoustic  Reverberaticm 
Special  Research  Project  (ARSRP).  The  study  area  is  shown  in  Figure  3;  it  occupies  75.000  kni^ 
on  the  west  flank  of  the  Mid-Atlantic  Ridge  n<»th  of  the  Kane  Fracture  Zone  and  comprises  crust 
with  ages  ranging  from  3  to  30  My.  The  underway  geqihysical  data  collected  in  digital  framat 
during  the  survey  include  Atlas  Hydrosweep  multibeam  bathymetry,  Hawaii  MRl  sidescan  sonar, 
magnetometer  and  gravimeter  profiles,  3.5-kHz  echosoundings,  and  single-channel  watergun 
profiles  [Tucholke  et  al.,  1993].  The  track  lines  were  run  at  an  average  speed  of  9.8  knots  on  a 
WSW-ENE  trending  grid  subparallel  to  plate  flow  lines  with  a  spacing  that  varied  from  4-6 
kilometers  near  the  ridge  crest  to  8-9  kilometers  in  deeper  water,  a  configuratimi  designed  to 
optimize  the  bathymetric  coverage  and  insonification  of  the  dominant  topographic  grain.  Global 
Posititming  System  navigation  was  available  throughout  the  survey. 

The  primary  data  enq)loyed  in  this  study  are  the  interpolated  depth  values  logged  each  minute 
from  the  center  beam  of  the  Hydrosweep  system.  The  average  along-track  spacing  between  data 
points  is  ~3(X)  m.  which  we  took  as  the  value  of  u  in  calculating  the  slope  statistics.  Based  cm  a 
detailed  examination  of  the  swath  maps,  the  track  segments  were  grouped  into  19  small  areas 
(1000-3000  km^)  where  the  statistical  character  of  the  abyssal  hills  was  more  or  less 
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homogeneous.  We  avoided  regions  where  the  abyssal  hills  wa«  interrupted  by  ridge-segment 
boundaries,  which  are  typified  by  deeper  bathymetry,  thicker  sediments,  and  inhomogeneous 
statistics  [Tucholke  et  al.,  1993].  The  track  lengths  in  each  area  varied  from  130  km  to  360  km, 
with  an  average  of  about  200  km.  Average  crustal  ages  were  determined  from  the  magnetic- 
anomaly  data  processed  by  M.  Tivey  [personal  conmunication,  1993]. 

//,  L,  and  ic  woe  estimated  for  each  area  by  fitting  the  center-beam  slope  distributions  using  a 
weighted  least-squares  algorithm;  the  results  plotted  as  numbers  on  Figure  3.  In  Figure  4,  we 
illustrate  how  the  slq)e-distribution  data  for  the  area  centered  at  26.0°  N,  46.5°  W  are  matched  by 
the  theoretical  model,  and  we  display  the  resulting  uncertainties  in  the  model  parameters  by 
contouring  their  95%  confidence  region.  The  best  estimates  are  //  -  220  ±7m,Ls34±3m,  and 
ir=  0.09  ±  .02  mVyr.  The  relative  sizes  of  the  standard  errors  (4%,  10%,  and  20%,  respectively) 
are  typical  of  the  otl^  areas.  The  confidence  re^cms  deviate  from  an  elliptical  shape,  reflecting  the 
non-Gaussian  character  of  the  estimation  uncertainties.  Their  basic  properties  can  be  easily 
understood  from  the  variations  shown  in  Figure  2.  The  strong,  positive  correlation  between  H  and 
L,  for  example,  is  related  to  the  smoothing  of  the  topography  by  the  sediments;  e.g.,  a  large  value 
of  5(0o4<)  may  correspond  either  to  thin  sediments  on  low-amplitude  basement  or  to  thick 
sediments  on  high-amplitude  basement.  There  is  essentially  no  tradeoff  between  L  and  k,  although 
the  95%  confidence  region  for  the  laner  is  distinctly  asymmetric.  The  upper  bound  on  the 
diffusivity  is  less  well  constrained  than  its  lower  bound,  because  the  part  of  the  slope  distribution 
sensitive  to  k  becomes  compressed  to  smaller  values  of  0  as  x  increases. 

With  these  uncertainties  in  mind,  we  turn  our  attention  to  the  regional  difrerences  evident  in 
Hgure  3.  H  varies  from  188  m  to  285  m,  with  a  mean  H  »  230  ±  7  m.  The  variation  is  not  large, 
but  it  does  correlate  with  the  ridge-segmoit  boundary  that  passes  through  the  center  ci  the  coiridor. 
The  average  abyssal-hill  amplitude  in  the  seven  areas  north  of  the  segment  boundary  is 
significantly  larg^  =  255  ±  12  m)  than  in  the  twelve  areas  to  the  south  {H  ^  214  ±  5  m), 
inqrlying  a  persistent  discontinuity  in  the  ridge-crest  processes  that  give  rise  to  the  abyssal-hill 
tc^giiq>hy.  In  contrast,  the  variation  in  L  is  primarily  in  an  east-west  direction,  along  the  plate 
flow  lines.  In  Figure  5,  we  plot  the  sediment  load  as  a  function  of  crustal  age.  It  is  thin  on  the 
young  crust  near  the  ridge  crest,  increasing  westward  to  maximum  thicknesses  as  high  as  55  m  on 
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crust  with  ages  of  10-15  Ma.  A  least-squares  regression  of  the  data  out  to  12  Ma  yields  an 
qjparent  sedimenoition  rate  of  3.4  nVMy. 

On  <dder  seafloor,  L  decreases  with  age,  droi^ing  to  values  as  low  as  20  m  at  25  Ma.  This 
observatioo  agrees  with  the  3  J-kHz  and  single-channel  seismic  data  collected  during  the  ARSRP 
survey.  A  detailed  discussion  of  the  trends  in  Figure  5  is  beytmd  the  scope  of  diis  initial  report, 
but  we  note  that  the  non-mcmotonic  character  of  L  can  plausibly  be  related  to  a  change  in  the 
carbcMuue  ctmqtensation  depth  (CCD).  Tucholke  and  Vogt  [1979]  have  hypotheazed  that,  prior  to 
die  late  Miocene,  the  (XD  in  the  central  North  Atlantic  was  shallower.  According  to  their  Hgure 
3,  die  most  rapid  change  in  the  CCD  occurred  in  die  mid-Miocene,  which  coincides  with  the 
sediment  peak  observed  in  our  data.  C^ubonates  dominate  the  upper-ocean  sediment  flux  in  this 
region,  and  their  rate  of  dissolution  would  have  been  greater  prior  to  this  time.  Taken  at  face 
value,  however,  the  nqnd  decrease  in  L  between  15  and  25  Ma  is  inconsistent  with  a  constant  (X!D 
and  sediment  flux  throughout  the  early  Miocene.  One  possibility  is  that  the  CCD  at  the  end  of  the 
Oligocene  was  even  shallower  than  Tucholke  and  Vogt  have  estimated,  and  that  the  period  over 
which  it  deq;)ened  included  the  eariy  Miocene. 

The  individual  areas  show  iqiparent  diffusivities  ranging  from  0.04  to  0.26  is^/yr.  Some 
variation  in  this  parameter  is  evidently  warranted  by  the  data;  attempts  to  fit  the  data  with  a 
cmistant  value  of  ir  led  to  unsatisfachny  results,  consistent  with  the  error  arutiysis  in  Hgure  4, 
which  indicates  that  the  observed  order-of-magnitude  change  should  be  resolvable.  The  physical 
significance  of  the  differences,  if  any,  is  unclear.  No  pattern  emerges  from  Figure  3,  for  example. 
The  mean  values  ncrth  and  south  of  the  segment  boundary  are  statistically  identical  (0.15  ±  .020 
and  0.12  ±  .019  m^^,  respectively),  as  are  those  east  and  west  of  AV  W  (0.13  ±  .017  and  0.13  ± 
.023  n^/yr).  The  estimates  of  L  and  ic,  like  these  of  L  and  /f,  are  unccnelated. 

More  inqrartant,  the  mean  value  of  the  apparent  diffusivity  obtained  from  the  cumulative  slope 
distributions  is  precisely  estimated  by  our  technique:  7  »  0.13  ±  .014  m^/yr.  To  our  knowledge, 
this  is  the  first  quantitative  estimate  of  sediment  diffusivity  in  a  deep-sea  environment  For 
comparison,  we  note  that  pelagic  sediments  in  the  central  North  Atlantic  are  transported  at  a  rate 
about  two  (nders  of  magnitude  greater  than  compacted  sediments  in  arid  regions,  whoo  K"  10~^ 
w^fyr  [e.g.  Hanks  et  aL,  1984],  but  several  orders  of  magnitude  less  rapidly  than  those  in  deltaic 
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and  near-shme  marine  settings,  which  yield  apparent  diffusivities  of  ~  10^  mVyr  [e.g.  Flemings 
and  Jtmlan,  1989]. 


DISCUSSION 

Our  methodology  obtains  good  estimates  of  the  sediment  load,  apparent  diffusivity,  and  rms 
amplitude  of  the  basemoit  topography  ftom  single-channel,  narrow-beam  echosounding  profiles  in 
rough  abyssal-hill  terrains  where  standard  sediment  profiling  techniques  cannot  be  applied. 
Although  the  diffusion  approximation  used  to  invert  for  these  parameters  undoubtedly 
oversimplifies  the  processes  that  govern  post-depositional  transport,  it  represents  a  first  step 
towards  a  dynamical  model  of  sedimentation  on  young  oceanic  crust.  Application  of  the  statistical 
techniques  presented  here  to  dense,  high-quality  data  sets  like  those  being  collected  in  the  ARSRP 
cofiidor  will  lay  the  foundation  for  future  improvements.  In  the  short  term,  however,  using  them 
to  process  existing  narrow-beam  data  may  be  helpful  in  mapping  the  distribution  of  sediments  on 
die  flanks  of  the  mid-ocean  ridges. 
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FIGURE  Captions 


Hg.  1.  Plan-view  relief  images  (top)  and  vertical  cross-sections  (bottom)  of  two  sedimented 
seafloors  obtained  by  the  diffusive-transport  algorithm  described  in  the  text,  showing  the  effect  of 
different  diffusivities;  ir*  0.02  m^/yr  (l^),  ir=  0.8  afifyr  (right).  The  models  have  the  same  age 
(7  B  25  Ma),  sedimentation  rate  (F  »  4  m/My).  sediment  load  (L  ~  100  m)  and  basement 
topography  (//  *  250  m,  X,  =  22  km,  =  5.4  km,  =  35®,  D  =  2.2).  The  lower  diffusivity 
produces  a  morphology  in  which  sediments  drape  hillsides,  while  the  higher  value  allows  a  mote 
rqrid  downslope  transport,  exposing  the  peaks  and  creating  flat  sediment  ponds.  The  dimensions 
of  the  images  are  30  km  x  30  km;  color  changes  at  elevations  of  -200, 200,  and  600  m.  Cross- 
sections  are  taken  along  the  dashed  lines  and  plotted  with  a  vertical  exaggeration  of  7.5 : 1.  The 
sectums  show  basement  topography,  the  sediment-water  interface,  and  two  intermediate  horizons 
correspcmding  to  isochrons  at  5  My  (L  =  20  m)  and  10  My  (50  m). 

Hg.  2.  Cumulative  slope  distribution  functions  S($,  u)  for  u  ^  300  m,  calculated  from  proHles 
aooss  sedimented-seafloor  models.  Top  diagram  illustrates  the  effect  of  varying  sediment  load  (L 
=  25, 50, 75, 100  m)  for  fixed  abyssal-hill  amplitude  and  diffusivity  (H  =  2(X)  m,  ir=  .10  m^/yr). 
Bottom  diagram  shows  the  effect  of  varying  diffusivity  (k=  .02,  .10,  .30,  .80  mVyt)  for  fixed  L 
(50  m)  and  H  (150  m  in  upper  curves,  325  m  in  lower  curves).  All  calculations  are  based  on  F  = 
4  m/My. 

Fig.  3.  Map  of  the  ARSRP  cotridor,  showing  bathymetry,  Ewing  9208  track  lines,  and  values  of 
die  model  parameters  estimated  for  the  19  areas  analyzed  in  this  report  Track  lines  used  to  obtain 
slope  distributions  are  highlighted;  L  (in  meters)  is  given  by  large  number  centered  on  each  area, 
widi  H  (in  merers)  and  k  (in  m^/yt)  listed  as  small  numbers  above  and  below.  Color  changes 
coirespcMid  to  water  dqiths  of  3200, 4000,  and  4800  m.  Region  of  deep  bathymetry  (light  brown) 
trending  WNW  through  the  center  of  the  ctnridor  is  a  ridge-segment  boundary  that  separates  hi^er 
an^litude  abyssal  hills  to  the  north  ( H  =  255  ±  12  m)  from  lower  amplitude  hills  to  the  south  ( H 
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■  214  ±  S  m).  The  sediment  load  L  ranges  from  1 1  to  55  m,  primarily  as  a  function  of  crustal  age 
(see  Figure  5);  maximum  values  occur  in  the  center  of  the  c(»rujk»^  (7  «  10-15  Ma).  Significant 
differences  in  the  apparent  diffusivity  are  obtained  fiom  the  data  analysis  (ir«  0.04-0.26  m^/yr), 
but  no  qratial  pattern  is  evident  Averaged  over  all  19  areas,  7  «  0.13  ±  .014  mVyr. 

Fig.  4.  7c{p  diagram  is  a  plot  of  5(6,  u)  for  the  area  centered  at  26.0°  N,  46.5°  W;  dots  are  the 
data  obtained  fiom  170  km  Hydrosweep  center-beam  souiulings,  and  the  curve  is  the  cumulative 
slope  distribution  doived  from  the  best-fitting  model  topograi^y(/f«  220  ±7  m,L«34±3tn,  k 
« .09  ±  .02  tr^/yr).  Lower  diagrams  are  L  vs.  i/  (center)  and  L  vs.  K(bottom),  showing  the  best 
estimates  (crosses)  and  95%  confidence  regions  (contour  lines). 

Fig.  5.  Estimates  of  sediment  load  L  as  a  function  crustal  age  for  the  19  areas  analyzed  in  this 
report  line  cortesptMids  to  a  constant  sedimentation  rate  of  3.4  m/My.  Shaded  region  delimits  the 
interval  of  CCD  deepening  discussed  by  Tucholke  and  Vogt  [1979].  We  speculate  that  the  non- 
monomiic  behavior  observed  in  the  ARSRP  corridor  is  due  to  this  change  in  the  CXD. 
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Modeling  Sedimentation  of  a  Rough  Seafloor 


Helen  F.  Webb 

Abstract  Sedimentation  is  the  primary  modifier  of  small-scale  tqx>graphy  created 
by  volcanic  and  tectonic  processes  in  the  vicinity  oi  ridge  crests.  In  particular,  sediment 
accumulation  tends  to  smooth  small-scale  topographic  roughness.  To  study  this  process  in 
pelagic  environments,  we  model  deep-sea  sedimentation  as  a  diffusion-like  process 
characterized  by  two  parameters,  a  sediment  flux  F  and  a  diffusivity  k.  We  did  not  account 
for  sediment  compaction  or  the  erosion  of  bedrock. 

For  low  diffusivity  values  (k<  20,000  m^/My)  our  model  causes  sediments  to 
drq>e  over  the  basement  topography.  For  high  values  (ir>  100,000  m^/My)  sediments  are 
swept  downslope  and  accumulate  in  topographic  lows  to  form  sediment  ponds. 

We  found  that  the  burial  of  a  hard-rock  surface  characterized  by  a  Gaussian  height 
distributitxi  generates  a  sediment-water  interface  whose  statistics  have  a  distinctly  non- 
Gaussian  character.  As  average  sediment  depth  increases,  the  RMS  height  and  fractal 
dimension  decrease,  and  the  characteristic  length  scales  of  the  topography  increase. 
Increasing  the  diffusivity  constant  of  the  sediment  distribution  process  also  decreases  the 
RMS  height  and  increases  die  characteristic  length  scale  of  the  topography. 


INTRODUCTION 

The  floor  of  the  world's  oceans  is  subject  to  a  constant  rain  of  sedimentation  from 
wind-carried  ccmtinental  dust,  the  weathering  of  oceanic  islands,  and  biological  sources. 
Althou^  sedimentation  rates  in  pelagic  environments  are  low  (from  less  than  1  m/My  to 
over  SO  m/My),  over  time  sedimentatitMi  snwoths  seafloor  topography.  Abyssal  plains, 
which  extend  oceanward  from  passive  continental  margins,  are  products  of  tens  of  millions 
of  years  of  sedimentatitm. 

The  processes  oi  sediment  transpon  and  deposition  in  the  abyssal  environment  are 
not  well  understood.  The  initial  location  of  sediment  deposition  must  be  a  function  of 
sediment  type  and  the  interaction  of  seabottom  currents  with  seafloor  constructs.  Over 
time,  steady  bottom  currents  [e.g.  Dickson  and  Kidd,  1987],  burrowing  animals,  and 
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localized  storms  (with  water  velocities  of  up  to  0.73  m/sec  [Richardstm  et  al.,  1981])  cause 
sediments  to  be  redistributed.  CSianges  in  txxtom  water  movement  due  to  global  climate 
changes  and  continental  reconfigurations  also  can  disturb  sediments,  leaving  gaps  in 
deepsea  drill  cores  [e.g.  Johnson,  1972;  Thiede  et  al,  1981].  These  processes  cause  the 
pattem  of  sedimentation  to  vary  between  different  locations.  In  some  regions,  sediments 
form  a  terrain-mimicking  blanket,  preserving  the  shape  of  abyssal  hills  created  at  the  ridge 
crest  (Figure  1).  In  other  regions,  sedimentatim  fills  up  valleys  while  leaving  the  tops  of 
abyssal  hills  bare,  quickly  evening  the  terrain  (Figure  2).  Most  areas  of  the  seafloor  fall 
between  these  extremes,  with  thicker  layers  of  sediment  covering  bedrock  lows  than  cover 
highs,  {holographs  of  the  seafloor  show  that  manganese  nodules  can  remain  uncovered, 
indicating  that  even  small  topographic  highs  can  be  kq>t  clear  of  sediments  [e.g.  Heezen 
and  Hdlister,  1971],  while  sediment  cores  show  that  cracks  between  pillow  basalts 
bectxne  filled  soon  after  creation  [e.g.  Reynolds  and  Linsley,  1985]. 

If,  over  time,  the  lateral  flux /of  erodible  sediment  is  proportional  to  the  local 
gradiait  of  height  via  a  difiusivity  constant  x;  sediment  transport  can  be  modeled  as  a 
diffusive  process  in  which  erodible  material  is  removed  from  topographic  highs  to  fill 
lows.  Cullings  [1960]  pioneered  this  approach  to  modelling  geomorphological  structures. 
Newman  and  Turcotte  [1990]  point  out  that  diffusion  will  produce  topography  that  is  not 
scale  invariant,  but  instead  preferentially  eliminates  short-wavelength  topography. 
Therefore,  they  conclude  that  diffusive  models  can  not  explain  most  of  the  topography 
found  in  nature.  Published  seafloor  power  spectra  fOT  ccmtinental  rise  regimes  tend  to  have 
in  increase,  instead  of  a  decrease,  in  power  at  high  wavenumber  [Fox  and  Hayes,  1985], 
so  it  is  unlikely  that  diffusion  would  adequately  explain  their  structure;  however,  the 
topognq>hy  of  the  abyssal  hill  regime  is  quite  different.  On  abyssal  hills,  the  main  effect  of 
sediment  movement  is  probably  the  filling  of  small  cracks,  which  would  reduce  power  at 
high  wavenumber.  Few,  if  any,  events  take  place  off  axis  over  time  which  would  increase 
power  at  high  wavenumber. 


Workers  have  successfully  used  diffusion  to  explain  the  evolution  of  scarplike 
landfonns  [Hanks  et  al.,  1984]  and  the  stratigraphy  found  in  foreland  basins  [Flemings  and 
Jordan,  1989]  and  in  shallow  marine  settings  [Kaufman  et  al.,  1991].  The  success  of  these 
models  is  due,  in  part,  to  the  incorporation  of  geological  processes  such  as  isostadc 
con^pensadm  [Flemings  and  Jordan,  1989]  and  the  decrease  of  wave  energy  with  depdi 
[Kaufman  et  aL,  1991].  In  our  model  of  sediment  transport,  we  have  incorporated  the 
geometry  of  the  seafloor  prior  to  sedimentation,  although  other  physical  factors 
contributing  to  sediment  distribution,  such  as  seabottom  currents,  have  been  ignored. 

The  total  height  h(x,t)  of  the  seafloor  is  given  by 

/i(x,r)  =  s(x,r)  +  b(x)  (1) 

where  six/)  is  the  thickness  of  the  erodible  layer  of  sediment  which  covers  a  durable 
basement  of  height  bix).  We  model  sediment  thickness  over  time  as  equal  to  the  gradient 
of/plus  the  downward  sediment  flux  F: 


dt  9x2  9x2 


subject  to  the  constraint  that  s(x,r)  ^  0.  We  compute  the  sediment  thickness  5(x,r)  by 


solving  (2)  numerically. 


SEDIMENTATION  ALGORITHM 


Numerical  methods  for  solving  parabolic  equations  are  based  on  the  approximation 

ds  _dh  _  hjXyt+At)  -  hix,t) 

df  (3) 

Values  of  s(x,0  are  calculated  by  iterating  forward  from  six,0)  =  0.  In  implicit  methods 

such  as  the  Crank-Nicolson  method  [Crank  and  Nicolson,  1947]  and  the  altemating- 

direction-implicit  scheme  [Peaceman  and  Rachford,  1955],  values  of  /i(x,r  -4-  dr)  are  not 

just  functions  of  /t(x,r)  but  also  involve  other  values  of  h  at  the  same  time  step.  These 
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methods  can  be  adapted  tt>  run  efficiently  (Xi  serial  con^uters.  However,  a  new  class  of 
computers,  including  the  Connection  Machine  manufactured  by  Thinking  Machines 
Coiporaticm,  is  built  to  perform  thousands  of  identical,  simple  operations  simultaneously 
on  grids  of  data.  Speed  is  increased  when  data  is  passed  only  between  neighboring  ^d 
points.  To  solve  (2)  using  such  a  machine,  we  use  an  explicit  method  so  that  values  of 
h(x4  +  At)aie  tmly  dependent  on  h(x,t\  allowing  all  values  of  h  for  each  time  step  to  be 
calculated  simultaneously.  Our  algorithm  is  a  variation  of  the  method  of  Schmidt  [1924], 
who  used  the  approximation 

_  r  ih(,x+Ax,t)  -  hjx^t))  +  ih(x-Ax,t)  -  h(x,r))l 


We  apply  (4)  and  (3)  to  (2),  obtaining,  fen:  the  two-dimensional  case, 

Jtoi+4()  s  -  /idrf))  +  (A(jr-ai,i)  -  *(W))]  +  differ) 

(6) 

subject  to  the  inequality  sOc^^At)  ^  0. 

The  speed  of  computation  is  proportional  to  the  magnitude  of  KAtUAxi^-,  however, 
for  xAtHAx^  >  1/3  for  the  two-dimensional  case,  and  icAtliAx)^  >  1/5  for  the  three- 
dimensional  case,  the  value  of  hix,t+At)  is  more  dependent  on  h  at  neighboring  points  than 
atx.  This  causes  h  to  oscillate.  For  1/3  <KAt/{Ax^  ^  1/2  for  the  two-dimensional  case, 
and  1/5  <KAtl{Ax)'^  ^  1/4  for  the  three-dimensional  case,  the  oscillations  damp  out  with 
increasing  t.  To  prevent  having  minimum  time  requirements,  and  to  maintain  conformity 
between  the  two-  and  three-dimensional  models  of  sediment  transport,  we  use 
jcdr/(ztc)2  =  0.2  for  both  algorithms.  Therefore,  the  timestep  At  needed  to  simulate  a 
given  diffusivity  k  is  calculated  by 

At  =  0.2(Ax)2/K.  (6) 
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For  a  given  number  of  years  7,  N  »  TIAt  iteradcms  are  needed 

Tbe  following  steps  comprise  the  two-dimensional  algorithm  This  algorithm  is 
easily  extended  to  three  dimensions.  A  graphical  example  of  the  diffusion  methodology  is 
provided  in  Figure  3: 

Initialization'.  Basement  topography  A(x)  is  discretized  on  a  straight  line  with 
horizontal  spacing  dx  and  total  horizontal  length  X  =  KAx. 

Execution  The  following  steps  are  then  performed  N  times: 

(a)  Sediment  Accumulation:  Parcels  of  dimensions  (AxJ^At)  are  accumulated 
at  each  horizontal  posititm  xk. 

(b)  Sediment  Redistribution:  First,  height  differences  Aht*  =  h(xii)  -  Kxk*i) 
and  Ahk'  =  h(xi^  -  hi^k-i)  are  calcula^  at  each  point.  To  conserve  sediment, 
d/to* »  0  and  d/tK'*’ «  0.  The  sediment  excess 

et  =  0.2(max(dAit*’,  0)  +  max(dAi*,  0))  (7) 

is  then  calculated.  If  d/tjfc+  >  0  and  sCxt)  <  ek  then  Ahk*  =  Ahk*'{sixk)feid  and 
Ahk+i’  is  set  equal  to  -1.0(d/ik''')-  Likewise,  if  dA^*  >  0  and  5(Xk)  <  ek  then 
Ahk'  =  Ahk'isixidlek)  and  dAk-i'*’  is  set  equal  to  -1.0(d/tk').  Finally,  the 
amount  of  sediment  change  for  each  point,  Ask  =  -0.2iAhk*‘  +  Ahk'),  is 
calculated  and  added  to  s(xid. 


MODEL  SEAFLOOR  BATHYMETRY 


We  use  Goff- Jordan  [1988]  model  seafloor  bathymetry  to  test  our  sedimentation 
methodology.  Model  bathymetry  was  used  instead  of  real  badiymetric  data  both  because 
we  lacked  sediment-free  data  with  uniform  point  spacings  in  the  25  - 100  m  range,  and  so 
that  we  could  more  easily  conq)aie  the  statistical  parameters  of  the  resulting  sedimented 
profiles  and  surfaces  to  the  original  tq)ography.  Goff  and  Jordan  [1988]  model  the 
seafloor  as  a  stationary,  zero-mean,  Gaussian  random  field  conqiletely  specified  by  its  two- 
point  covariance  function.  This  function  has  five  parameters  which  crarespond  to  RMS 
height  //,  characteristic  along-  and  across-flowline  comer  wavenumbers  and  ks ,  local 
strike  and  fractal  dimension  D. 

EXPERIMENTAL  DESIGN 

We  modelled  the  sedimentation  of  both  a  three-dimension  grid  of  simulated 
bathymetry  and  five  two-dimension  profiles  of  simulated  bathymetry.  Although  three- 
dimensional  modelling  probably  produces  the  more  realistic  results,  large  grids  are  required 
to  cover  multiple  abyssal  hills  using  even  a  course  point  spacing.  Reducing  Ac  by  a  facmr 
of  a  increases  both  the  grid  size  and  the  number  of  iteraticms  N  by  factors  of  a^,  leading  to 
a  total  computation  time  increase  of  oc^,  whereas  a  similar  decrease  in  point  spacing  along  a 
two-dimensional  profile  only  leads  to  a  total  confutation  time  increase  of  a^.  Therefore, 
two-dimensional  profiles  were  used  to  test  the  sedimentation  model  on  more  closely  spaced 
simulated  bathymetry. 

The  three-dimensional  grid  of  simulated  bathymetry  (Figure  4)  was  created  by 
tix>lding  a  surface  generated  using  the  Goff-Jordan  methodology  with  DBDB5  data  fiom 
the  vicinity  of  22®  N,  47®  W  (near  the  Mid  Atlantic  Ridge,  south  of  the  Kane  Fracture 
Zone)  [Goff  and  Jordan,  1990].  This  model  seafloor  bathymetry  is  already  available  to  the 


acoustics  community  via  the  Office  of  Naval  Research's  Special  Research  Project  on 
Acoustic  Reverberation.  It  consists  of  1000  x  lOUO  points  of  bathymetry  with  a  grid 
spacing  Ax  of  100  m.  This  spacing  was  chosen  so  that  several  abyssal  hill  wavelengths 
would  be  included  in  all  dircctions.  The  region  was  sedimented  for  20  My  of  model  time 
using  diffusivides  of  1600  m^/My,  16,000  m^/My,  160,000  m^/My,  and 
1,600,000  m^/My.  Flux  was  kept  constant  at  10  m/Myr,  and  sediment-water  interface 
surfaces  were  recorded  every  S  My  of  model  dme.  Several  of  the  surfaces  generated  are 
shown  in  Figure  5,  and  profiles  taken  across  the  middle  of  these  regicms  are  shown  in 
Figure  6. 

To  test  our  model  on  a  more  closely  spaced  data  set,  we  modeled  the  sedimentation 
of  five  125  km-long  profiles  which  had  relatively  small  horizontal  point  spacings  Ax  of 
25  m.  The  profiles  were  each  sedimented  for  10  My  of  model  time  using  diffusivides  of 
1600  m^/My,  16,000  m^/My,  160,000  m^/My,  and  1,600,000  m^/My.  Flux  was  kept 
constant  at  10  m/Myr,  and  sediment-water  interface  profiles  were  recorded  every  2.5  My. 
These  profiles  are  shown  in  Figure  7. 

The  parameters  governing  the  diffusive  process,  v,  total  time  T,  and  flux  F,  can  be 
combined  to  form  a  non-dimensional  measure  of  diffusivity 

TF^  (8) 

Thus,  for  a  given  average  sediment  thickness  (TF),  increasing  F  and  v  by  the  same  factor 
produces  equivalent  topography. 

RESULTS  AND  ANALYSIS 

We  found  that  our  model  was  able  to  produce  sediment  distributions  that  resemble 
those  found  in  nature.  Using  a  diffusivity  jcof  1.6  x  10^  m^/My  (Figure  5(a,b); 

Bgure  6(a);  Rgure  7(a)),  we  produced  topography  that  resembles  the  "pelagic  drape" 
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character  (tf  Figure  1.  Higher  values  of  jclead  to  increased  mobility  of  sediments,  with 
sediments  ponding  between,  and  eventually  covering,  abyssal  hills  (Figure  5(c,d); 

Figure  6(d):  Figure  7(d)).  This  topography  resembles  that  of  Figure  2. 

When  irand  t  are  small,  sediments  to  not  spread  far  from  where  they  are  originally 
located,  so  die  variation  of  sediment  thickness  has  a  low  amplitude  and  a  small 
characteristic  horizontal  scale  (Hgure  8(a)).  As  T  increases,  this  characteristic  scale 
increases  (Figure  8(b)),  as  sediments  are  gradually  able  to  fill  topographic  lows. 

Increasing  ir (Figure  8(c))  causes  both  the  arrqilitude  and  characteristic  scale  to  increase, 
owing  to  die  formation  of  large,  deep  sediment  ponds.  When  the  sediment  load  is  large 
enough  (Hgure  8(d)),  the  variations  in  r (x,/)  mimic  those  of  the  basement  topography  h(x) 
(Figure  4). 

The  effects  of  this  diffusivity-govemed  spread  of  sediment  can  be  seen  in 
histograms  of  sediment  thickness  and  seaflocx*  height  In  a  regime  of  high  sediment 
diffiisivity,  die  seafloor  height  distribution  quickly  loses  its  Gaussian  character  (Hgure  9). 
For  large  sediment  loads  (T  =  20  My),  the  height  distribution  contains  many  peaks, 
ctmre^xxiding  to  the  damming  of  sediments  bdiind  abyssal  hills.  Low-lying  areas  are 
quickly  filled,  causing  an  increase  in  the  minimum  heights  with  7;  however,  since  peaks 
do  not  accumulate  sediment  the  upper  tails  of  all  of  the  distributions  are  similar. 
Ifistograms  showing  the  height  distribution  of  sedimented  seafloor  for  7  =  20  My  and  a 
variety  of  diffiisivities  show  the  development  of  this  peakiness  (Figure  10). 

Sediment  thickness  histograms  (Hgures  1 1  and  12)  show  that,  where  sediment 
exists,  it  exhibits  a  Gaussian-like  distribution  of  depths.  However,  as  diffusivity 
increases,  much  of  the  seafloor  is  stripped  of  its  blanket  of  sediment 

We  used  the  results  of  sedimenting  long,  densely  spaced  (25  m)  two-dimensional 
profiles  shown  in  Hgure  7  to  determine  the  effects  of  sedimentation  on  seafloor  power 
spectra  and  the  Goff- Jordan  statistical  parameters.  The  power  spectra  of  sedimented 
profiles  resemble  the  power  spectrum  of  the  basement  topography  at  low  wavenumbers. 


but  feature  roll-offs  at  intennediate  wavenumbers.  At  high  wave  numbers,  the  slope  of  the 
power  spectra  of  the  sedimented  profiles  is  usually  the  same  as  the  slope  of  the  basement's 
power  spectrum. 


After  2.5  My  of  sedimentation,  the  greatest  deviations  from  the  spectrum  of  the 
Goff- Jordan  basement  tqwgraphy  occur  for  low  values  of  k,  where  the  sediment 
distribution  is  in  the  form  of  pelagic  drape  (Hgure  13).  The  smoothing  this  entails  causes 
low  power  at  high  wavenumbers.  At  high  values  of  k,  uncovered  peaks  contribute  energy 
at  high  wavenumbers.  After  10  My  of  sedimoitation,  all  diffusivities  yield  the  same  power 
at  high  wavenumber;  however,  increasing  irreduces  power  at  intermediate  wavenumbers 
(Figure  14). 

Plots  of  power  spectra  of  profiles  with  similar  diffusivities  show  that,  as  the 
amount  of  sediment  increases,  power  at  intermediate  and  high  wavenumbers  decreases 
(Figures  IS  to  18).  The  pattern  is  similar  fOT  all  diffusivities,  with  the  primary  difference 
being  a  downward  progression  of  the  wavenumber,  with  increasing  diffiisivity,  at  which 
the  power  spectra  of  the  sediment  profiles  separate  from  the  power  spectrum  of  the 
basement  topography. 

The  power  spectra  roll-offs  seen  are  ccnnmon  to  models  of  mass  distributicm  that 
entail  diffiisivity  [Newman  and  Turcotte,  1990].  Just  as  there  is  a  low  cmier  wavenumber 
in  the  power  spectrum  of  the  simulated  basement  topography  that  correspcmds  to  a 
characteristic  upper  topographic  wavelength,  the  bends  in  the  power  spectra  at  intermediate 
wavenumber  may  correspond  to  a  characteristic  wavelength  of  topography  that  has  been 
buried  ot  filled  in  by  sediment 

Although  the  power  spectra  of  the  sedimented  profiles  clearly  do  not  follow  the 
Goff-J(»dan  model,  fitting  ideal  covariance  functions  to  the  covariance  functions  of  the  data 
to  obtain  values  for  RMS  //,  characteristic  wavenumber  ,  and  fractal  dimensicm  D 
allows  us  to  quantify  the  effects  of  sedimentation  on  the  larger  scale  statistics  of  the 
basement  topography.  These  parameters  are  listed  in  Table  1.  The  immediate  decrease  in 
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finctal  dimensioo  shows  that  much  of  the  small-scale  characteristic  roughness  of  the 
topography  is  filled  in  by  sediment  soon  after  abyssal  hill  fonnation;  the  increase  in  D  with 
high  diffiisivity  is  caused  by  the  reemergence  of  peaks.  The  decrease  in  H  with  sediment 
thickness  (and,  dierefore,  time)  is  caused  by  valleys  filling  with  sediment,  reducing  abyssal 
hiU  height  Higher  difiusivities  probably  cause  sediment  particles  to  be  more  rapidly 
transported  from  peaks  to  values,  lowering  H  further.  The  decrease  in  characteristic 
wavenumber  kn  with  sediment  thickness  and  diffiisivity  corresponds  to  an  increase  in  the 
characteristic  length  scale  of  the  topogrt^hy,  which  is  primarily  caused  by  the  submergence 
of  lower-lying  abyssal  hills. 


CONCLUSIONS 

We  have  developed  a  model  for  the  sedimentation  of  abyssal  hills.  Sediment 
transport  is  modeled  as  a  diffusive  process  in  which  mass  nxives  relative  to  local  gradients. 
I¥ofiles  and  surfaces  created  using  our  methodology  vary  with  the  diffiisivity  constant,  but 
have  the  same  diaracter  as  sections  of  the  actual  seafloor  (Bguies  1  and  2). 

Our  model  changes  the  depth  distribution,  power  spectra,  and  statistical  parameters 
of  die  seafloOT.  In  particular,  characteristic  abyssal  hill  heights  are  decreased  while  the 
distances  between  peaks  are  increased  A  comer  in  power  spectra  at  intermediate 
wavenumber  may  indicate  an  lower  limit  for  the  wavelength  of  topographic  features. 

Our  model  is  undoubtedly  an  oversimplificaticMi  of  the  processes  that  govern 
sedimentation.  For  instance,  we  assiune  that  there  is  no  preferential  directicm  for  the  lateral 
movement  of  sediment  particles;  however,  in  nature,  unidirectional  seafloor  currents 
probably  do  affect  sediment  placement  We  also  assume  that  diffusivity  is  constant 
throughout  a  region;  however,  sedimentation  processes  may  differ  between  peaks  and 
valleys.  To  investigate  how  the  sedimentaticxi  of  an  actual  seafloor  differs  from  our  model, 
we  plan  to  apply  our  methodology  to  grids  of  estimated  basement  depths  and  compare  the 


resulting  sedunented  surface  to  actual  bathymetric  data.  Detailed  data  to  be  used  for  this 
comparison  is  currently  being  collected  as  part  of  a  survey  of  the  western  flank  of  the  Mid 
Atlantic  Ridge  south  of  the  Kane  Fracture  Zcme.  This  data  set  will  also  be  used  to  develop 
methods  for  invening  single-channel  bathymetry  to  determine  the  sediment  diffusivity 
parameter,  which  may  be  a  useful  descriptor  of  the  style  of  sedimentation  of  the  seafloor. 
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HGURE  CAPTIONS 


Figure  1:  These  seiignic  reflection  profiles  fiom  the  North  Atlantic  southwest  of  the  Kane 
Racture  zone  show  sediments  draping  topography.  Estimated  diffusivities  for  the  top  two 
profiles  are  <  10,000  m^/My.  The  bottom  profile  has  an  estimated  diffusivity  of 
50,000  m^/My.  Verticle  exaggeration  is  approximately  10:1. 

Hgure  2:  The  sediments  shown  in  this  seismic  pix^e  have  formed  flat  sediment  ponds. 
The  estimated  diffusivity  for  this  regitm  is  1,000,000  m^/My.  This  profile  images  crust  on 
die  eastern  side  of  the  Mid  Atlantic  Ridge  soudi  of  the  Kane  Fracture  Zjoot,  Verticle 
exaggeration  is  approximately  10:1. 

Figure  3:  A  sample  basonent  topography  is  subjected  to  two  iteratitxis  of  the  sediment 
distribution  algorithm.  The  quantity  FAx  here  is  equal  to  O.S  m,  and  K-3.  First, 
basement  topography  is  discretized,  yielding  b(0)  =  1.0,  b(l)  =  3.0,  bi2)  =  2.0,  and 
6(3)  *  2.0.  At  f «  0,  h(x)  *  b(x)  (box  (a)). 

After  adding  FAx  to  each  horizontal  position  xjt;,  6(0)  =  1.5, 6(1)  =  3.5, 6(2)  =  2.5,  and 
6(3)  =  2.5  (box  (b)). 

Following  the  procedure  described  in  the  text,  height  differences  are  then  calculated,  giving 
Aho-  =  0.0,  Aho*  =  -2.0,  A6r  =  2.0,  A6i+  =  1.0,  A62*  =  -1.0,  A62+  =  0.0,  A63- 
=  0.0,  and  A63+  =  0.0.  From  equation  (8),  we  then  calculate  eo  =  0.2(0.0)  =  0.0,  ei  = 
0.2(3.0)  *  0.6,  €2  =  0.2(0.0)  =  0.0,  and  €3  -  0.2(0.0)  =  0.0.  Since  Ahi*  >  0  and  j(xi) 
<  Cl  a  new  Ahi*  =  d6i+(j(xi)/ci)  is  calculated,  giving  Ahi*  =  1.0(0.5/0.6)  =  0.833. 
462*  is  then  reset  to  -0.833.  Also,  since  Ahi'  >  0  and  s(xi)  <  ei  a  new  46r  =  zl6r 
(j(xi)/ci)  is  calculated,  giving  46r  =  2.(K0.5/0.6)  =  1.667.  Aho*  is  then  reset  to  -1.667. 
Using  Ask  =  -0.2(46*+  +  46*').  we  then  compute  Asq  =  0.333,  Asi  =  -0.5,Aso  = 

0.167,  and  Aso  =  0.0.  Adding  these  changes  in  sediment  amount  to  s(Xk)  yields  6(0)  = 
1.833,  6(1)  =  3.0,  6(2)  =  2.667,  and  6(3)  =  2.5  (box  (c)). 


This  completes  (me  iteration  of  the  algorithm.  Repeating  the  steps  above,  adding  FAx  to 
each  horizontal  positicm  produces  the  configuration  in  box  (d).  After  redistribution,  = 
2.567,  A(l)  =  3.2,  h{2)  =  3.2,  and  /i(3)  =  3.033  (box  (e)). 

Figure  4:  Grey-shaded  color  relief  image  of  the  basement  topogrt^hy  b(x)  employed  in  our 
sediment  accumulation  calculations.  This  model  was  created  by  molding  a  surface 
gmerated  using  the  Goff- Jordan  [1988]  methcxiology  with  DBDBS  data  frcxn  the  vicinity 
of  22°  N,  47°  W  (near  the  Mid- Atlantic  Ridge,  south  of  the  Kane  Bacture  Zone).  Grid 
spacing  is  100  m.  Goff-Jordan  parameters  used  to  create  the  interpolating  surface  are 
H  =  225  m,D  -  2.2,  kn  =  0.55  km~^  (corresponding  to  a  characteristic  abyssal  hill 
spacing  of  5.9  km)  and  kj  =  0. 15  km'^  (corresponding  to  a  characteristic  abyssal  hill  length 
of  22  km). 

Figure  5:  Maps  of  total  topography  of  the  seafloor  h(x,t)  after  sedimentation,  showing 
four  cases: 


(a) 

r  =  5  My, 

ic=  1.6  X  10^  m^/My 

(b) 

f  =  20  My, 

*•=  1.6  X  lO^m^/My 

(c) 

f  =  5  My, 

K=  1.6  X  10^  m^/My 

(d) 

r  =  20My, 

v=1.6xl0^m^/My 

In  all  four  cases,  the  sediment  flux  F{x,t)  was  held  constant  at  10  m/My. 

The  low-diffusivity  cases  (a)  and  (b)  result  in  "pelagic  drape",  with  an  essentially  uniform 
thickness  of  sediments  overlying  the  basement  tc^graphy;  h{x,t)  therefore  resembles  a 
slightly  sm<x>thed  version  of  b(x). 


For  die  high-diffiisivity  cases  (c)  and  (d),  die  increased  mobility  of  the  sediments  foims 
sedinoent  ptmds  between  abyssal  hills;  after  20  My  of  sediment  accumulation  (d),  only  the 
Molest  peaks  emerge  from  the  sediment 

Figure  6:  Cross-sections  through  three-dimensional  sedimented  seafloor  nKxlels,  showing 
the  variations  in  iscqiachs  due  to  variations  in  the  difriisivity  xi  The  black  lines  represent 
the  basement  topography  corresponding  to  an  east-west  profile  through  the  center  of 
Figure  4.  The  colored  lines  are  isopachs  with  spacings  of  S  My.  Diffusivities  used  to 
compute  the  profiles  are: 

(a)  JC=  1.6  X  10^  m^/My 

(b)  »r=  1.6  X  10^  m^/My 

(c)  fc=  1.6  X  10^  m^/My 

(d)  ic=1.6xl0^m^/My 

Venical  exaggeration  is  10:1. 

Hgure  7:  Portitms  of  model  basement  (black)  and  computed  sediment  isopachs  (red)  from 
one  of  the  densely  spaced  (25  m)  two-dimensional  profiles  used  for  determining  the  effects 
of  sedimentation  on  seafloor  power  spectra  and  Goff- Jordan  parameters.  The  basement 
profiles  have  Goff- Jordan  parameters  H  =  1(X)  m,D  =  1.4,  and  k„  =  0.4  km'^ 
(correqxmding  to  a  characteristic  abyssal  hill  spacing  of  7.4  km).  Isopach  spacing  is 
2 J  My.  Diffusivities  used  to  compute  sediment  profiles  are: 


Vertical  exaggeradon  is  10:1.  These  diagrams  are  scaled  so  that  abyssal  hill  heists  and 
sediment  layer  thicknesses  are  approximately  the  same  as  in  Hgure  4. 

Figure  8:  Maps  of  sediment  thickness  s(x,t),  corresponding  to  the  four  cases  (a)-(d)  in 
Figure  S.  For  all  of  these  plots,  die  salmon-red  division  of  die  color  scale  occurs  at  the 
avenge  sediment  thickness,  which  is  SO  m  for  r  =  S  My  (a,  c)  and  200  m  for  r  =  20  My 
(b,  d). 

When  xrand  t  are  small  (case  a),  the  variation  of  sediment  thickness  has  a  low  amplitude 
and  a  small  characteristic  horizontal  scale.  As  r  increases,  this  characteristic  scale  increases 
(case  b).  Increasing  x'(case  c)  causes  both  d^  an^limde  aixl  characteristic  scale  to 
increase,  owing  to  the  fonnation  of  large,  deep  sediment  ponds.  When  the  sediment  load  is 
large  enough,  the  variations  in  5(x,r)  mimic  those  of  the  basement  topography  b(x). 

Hgure  9:  Histograms  showing  the  hei^t  distributimi  of  the  basement  topography  (blade 
line)  and  sedimented  seafloor  for  r  =  5  My  (red),  10  My  (blue),  15  My  (green),  and  20  My 
(purple),  confuted  for  it  =  1.6  x  10^  m^/My  (high  diffusivity).  For  large  sediment  loads 
(r  =  20  My),  the  height  distribution  has  many  sharp  peaks.  Low-lying  areas  are  quickly 
filled,  causing  an  increase  in  the  minimum  heights  with  r,  however,  since  peaks  do  not 
accumulate  sediment,  the  uppo*  tails  of  all  of  the  distributions  are  similar. 

Hgure  10:  Histograms  showing  the  height  distribution  of  sedimented  seafloor  for 
r  s  20  My  and  *:=  1.6  x  10^  m^/My  (red),  1.6  x  10^  m^/My  (blue),  1.6  x  10^  m^/My 
(green),  and  1.6  x  10^  m^/My  (purple).  Note  the  development  of  the  peaks  discussed  in 
Hgure  9. 

Hgure  11:  Histograms  showing  the  distribution  of  sediment  thicknesses  for  r  =  5  My 
(red),  10  My  (blue),  15  My  (green),  and  20  My  (purple),  computed  for  )c=  1.6  x  10^ 
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m^/My  (high  difiusivity).  The  percentage  of  seafloor  covered  by  less  than  one  meter  of 
sediment  is  39%,  25%,  17%,  and  1 1%,  resp^vely. 

Figure  12:  I£stograms  showing  the  distributitm  of  sediment  thicknesses  fm:  r  s  20  My 
and  If  =  1.6  X  10^  m^/My  (red),  1.6  x  10^  m^/My  (blue),  1.6  x  10^  m^/My  (green), 
and  1.6  x  10^  m?/My  (purple).  The  percentage  of  seafloor  covered  by  less  than  one 
meter  of  sediment  is  0%,  0.0003%,  3.5%,  and  12%,  respectively. 

Hgure  13:  The  power  spectra  of  the  2.5  My  layer  subject  to  diffiisivities  of 
1.6  X  10^  m^/My,  1.6  x  10^  m^/My,  1.6  x  10^  m^/My,  and  1.6  x  10^  m^/My. 

Hgure  14:  The  power  spectra  of  the  10  My  layer  subject  to  diffusivities  of 
1.6  X  10^  m^/My,  1.6  x  10^  m^/My,  1.6  x  10^  m^/My,  and  1.6  x  10^  m^/My. 

Hgure  15:  The  power  spectra  of  the  basement  and  layers  recorded  after  2.5, 5, 7.5,  and  10 

a  7 

My  of  sedimentation  subject  to  a  diffiisivity  of  1.6  x  10“^  m'^/My. 

Figure  16:  The  power  spectra  of  the  basement  and  layers  recorded  after  2.5, 5, 7.5,  and  10 
My  of  sedimentation  subject  to  a  diffiisivity  of  1.6  x  10^  m^/My. 

Figure  17:  The  power  spectra  of  the  basement  and  layers  recorded  after  2.5, 5, 7.5,  and  10 
My  of  sedimentation  subject  to  a  diffiisivity  of  1.6  x  10^  m^/My. 

Hgure  18:  The  power  spectra  of  the  basement  and  layers  recorded  after  2.5, 5, 7.5,  and  10 
My  of  sedimentation  subject  to  a  diffiisivity  of  1.6  x  10^  m^/My. 


Table  1 


)f(m2/My) 

r(My) 

tfi  (m^) 

kt,  (cyclcs/km) 

Fractal  E>imei 

0 

0 

9538 

0.45 

1.4 

1600 

2.5 

9326 

0.55 

1.2 

1600 

5.0 

9212 

0.6 

1.1 

1600 

7.5 

9124 

0.6 

1.05 

1600 

10.0 

9050 

0.6 

1.05 

16,000 

2.5 

8923 

0.6 

1.05 

16,000 

5.0 

8506 

0.55 

<1.05 

16,000 

7.5 

8210 

0.55 

<1.05 

16,000 

10.0 

7976 

0.5 

<1.05 

160,000 

2.5 

8389 

0.55 

1.05 

160,000 

5.0 

7323 

0.45 

<1.05 

160,000 

7.5 

6531 

0.4 

<1.05 

160,000 

10.0 

5947 

0.35 

<1.05 

1,600,000 

2.5 

8053 

0.5 

1.05 

1,600,000 

5.0 

6458 

0.35 

1.1 

1,600,000 

7.5 

5093 

0.3 

1.1 

1,600,000 

10.0 

4054 

0.25 

1.1 

a 


b 


kUomtna 


Basement 


UUmtUn 


height  (m) 


Figure  5(t) 


Seafloor,  K=1600m2/My  T^20My 


m 


-sooi)  -4000  -:ioao  -0000 


Figure  5(c) 


Seafloor y  k  =  1,600,000  m^/My  7  =  5  My 


ZO  40  60  80  100 


kUomeurs 


~‘n){):)  -fOOG  SQGD 


height  (m) 


kUsmtUn 


Seafloor,  K=  1,600,000  n^lMy  T  =  20My 


height  (m) 


10G 


hBometers 


O 

O 


o 

o 


(f>)^ 


>N 


£ 


o 

o 

o 

o 

o 

CO 

0k 


II 

Vci 


sja:^8LU 


i 


i 


i 


i 


i 


i 


( 


< 


( 


i 


(y)l 


f 


sjapuu 


I 


2D  Profile:  k  =  16,000  mV^y 


SJ9|9UU 


i 


kUeauttn 


Sediment  Thickness,  K-ldOOm^lMy  T-SMy 


kilometers 


thickness  (m) 


kUoHUtm 


I 


Sediment  Thickness,  1600  m^lMy  T- 20  My 


thickness  (m) 


Fi'0t/re  6^c) 


Sediment  Thickness,  1,600,000  m2 /My  T  =  5  My 


40  00  SO  100 


kOometen 


thickness  (m) 


Seafloor  Height:  k  =  1,600,000  m^/Myr 


4400  -4000  -3600  -3200  -2800 

height,  meters 


Seafloor  Heig 


OI 


S0ou9jjnooo 


ight,  meters 


Sediment  Thickness:  k  =  1,600,000  m^/My 


saouajjnooo 


200  300  400  500  600 

thickness,  meters 


Sediment  Thickness,  T  =  20  My 


S80U9Jjn000 


100  200  300  400 

thickness,  meters 


2D  Algorithm 


WAVENUMBER,  cycles/km 


2D  Algorithm:  T  =  10  My 


-hi 


WAVENUMBER,  cycles/km 


Algorithm:  k  =  1 ,600  roP-/Uy 


WAVENUMBER,  cycles/km 


2D  Algorithm:  /c  =  16,000  roP-/Uy 


(uj>i/s8ioAo)/2Ui  •y3M0d 


WAVENUMBER,  cycles/km 


Algorithm:  k  —  1,600,000  m^/My 


WAVENUMBER,  cycles/km 


Stochastic  interpolation  of  a  Wide-Beam  Bathymetric 
Profile  Within  the  ONR-ARSW^  Acoustic 
Reverberation  Corridor 


John  A.  Goffl,  Helen  F.  Webb^,  and  Thomas  H.  Jordan^ 


^Department  cf  Geology  and  Geophysics. 

Woods  Hole  Oceanographic  Institution,  Woods  Hole.  MA  02543 

^Department  of  Earth.  Atmospheric  and  Planetary  Sctnces. 
Massachusetts  Institute  of  Technology,  Cambridge,  MA  02139 


SUMMARY 

This  technical  report  describes  the  procedure  implemented  in  the  stochastic  interptdation 
of  a  bathymetric  profile  generated  by  a  wide-beam  echosounder.  The  principal  purpose  of 
this  interpolation  is  to  provide  acoustic  modelers  with  realistic,  densely  sampled  profiles 
from  which  to  predict  backscatter  phenomena.  The  profile  chosen  for  interpolation, 
collected  on  the  Verna  2602  cruise  (V2602),  passes  thrwgh  the  ONR-ARSRP's  acoustic 
reverberation  corridor  in  the  North  Atlantic.  After  gridding  the  entire  profile  at  data 
spacing  of  1  km,  three  20  km  subsets  were  chosen  for  stochastic  interpolation  to  1  m 
data  spacing.  These  sites  are  located  near  target  areas  for  the  1991  acoustic  reconnaissance 
mqierimenL  Two  sites,  B  and  C.  cotieqxmd,  respectively,  to  a  rough  plateau  and  a  steep 
escarpment  Site  A  corresponds  to  a  sediment  pond.  For  this  site  we  estimated  basement 
dqiths  from  the  V2602  seismic  profile  and  applied  the  interpolation  to  basement  depths 
only.  For  sites  B  and  C  we  postulate  a  thin  veneer  of  sediments  which  fill  the  sm^ler 
nooks  and  crannies  and  which  cannot  be  detected  by  the  seismic  inofiler.  To  simulate  this 
possible  effect  we  have  employed  a  numerical  model  for  sedimenting  a  rough  surface. 
The  original  and  interpolated  profiles,  along  with  sediment  depths  (observed  for  site  A. 
numerically  estimated  for  sites  B  and  C)  for  each  will  be  provided  to  all  interested 
ARSRP  participants. 
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1.  Introduction 


One  of  the  primary  goals  of  ONR's  Acoustic  Reverberation  SRP  (ARSRP)  is  to  model 
the  reveiberation  of  acoustic  energy  from  realistic  representations  of  the  ocean  floor.  The 
ocean  floor,  especially  the  basement  generated  at  mid-ocean  ridges,  is  rough,  and  contains 
features  at  all  measurable  scales.  Currently,  our  principal  and  most  successful  means  of 
modeling  acousdc  interaction  with  a  rough  intoface  is  through  numerical  modeling  (e.g. 
finite  differences,  finite  elements,  integral  equations,  etc.).  This  presents  us  with  an 
unfortunate  problem:  our  ability  to  deterministically  map  the  seafloor  falls  far  short  of  the 
detail  that  acoustic  modelers  need  to  accurately  predict  acoustic  reverberation.  For 
exaiiq)le,  for  the  wavelengths  considered  in  the  ARSRP,  -15  m,  the  seafloor  must  be 
accurately  represented  at  the  meter  scale  or  less.  By  contrast,  surface  operated  swath 
mapping  systems  (e.g.  Sea  Beam,  Sea  Beam  2000,  Hydrosweep,  SeaMARC  II,  etc.),  the 
only  bathymetric  instruments  that  have  the  ability  to  ensonify  large  enough  tracts  of  ocean 
floor  to  be  of  interest  to  Navy  needs,  have  lateral  resolutions  in  deep  water  of  100  m  at 
best 

For  some  purposes,  the  two  orders  of  magnitude  gap  between  what  is  needed  by 
acoustic  modelers  and  what  can  be  delivered  by  seafloor  mappers  can  be  filled  by 
extrapolation  to  high  wavenumbers  (interpolation  to  small  scales)  through  the  application  of 
stochastic  modeling  techniques.  Stochastic  models,  such  as  parameterizations  of  the 
covariance  function  or  power  spectrum,  even  very  simple  ones,  inherently  describe  the 
statistical  properties  of  morphology  at  all  scales.  For  example,  fractal  models  [e.g.  Bell, 
1975;  Mandelbrot,  1983;  Fox  and  Hayes',  1985;  Gilbert  and  Malinverno;  1988;  Goff  and 
Jordan;  1988]  predict  that  the  statistical  characteristics  of  small  scale  features  are  related,  ad 
infinitum,  to  those  of  large  scale  features  by  a  simple  self-affine  relationship.  Thus,  if  the 
large  scale/low  wavenumber  morphology  is  known,  the  small  scale/high  wavenumber 
moqriiology  can  be  predicted.  In  the  wavenumber  domain,  simple  fractal  models  can  be 
represented  as  straight  lines  in  a  log-log  plot  of  the  power  spectrum  for  a  topographic 
profile.  The  slope  of  this  line  determines  the  fractal  dimension.  High  wavenumber 
extrtqmlation  is  accomplished  simply  1^  extending  that  line  to  higher  wavenumbers  than  are 
resolved. 

In  the  abyssal  hill  morphology  encompassed  by  the  ARSRP  acoustic  reverberation 
corridor  (Figure  1),  we  have  little  indication  as  to  whether  the  statistical  behavior  of  small 
scale  features  can  accurately  be  predicted,  either  by  fractal  or  other  models,  from  the 
statistical  behavior  of  larger  scale  features.  The  geophysical  reconnaissance  portion  of  die 
ARSRP  will  address  this  question.  Some  encouragemoit  ccnnes  from  existing  data  in  other 
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r^ons,  however.  For  example,  using  both  Sea  Beam  and  Deep-Tow  bathymetric  data 
near  the  axis  of  the  Galapagos  ridge,  Goff  and  Kleinrock  [1991]  demonstrated  that  a 
unifonn  fractal  bduvior  can  extend  from  the  Idlometer  scale  down  to  the  ten  meter  scale. 

Realistic  stochastic  models  of  seafloor  morphology  aid  acoustic  modelers  in  two  ways; 
(1)  by  providing  the  basis  for  analytic  acoustic  models  on  ensembled  sets  and  (2)  by 
providing  the  means  of  generating  realistic  simulations  or  interpolations  of  batiiymetry  to 
arbitrary  resolution  and  scale  for  use  in  numerical  modeling  procedures.  It  is  in  the  latter 
capacity  that  tiiis  technical  report  is  directed.  In  the  first  of  these  ARSRP  technical  reports 
[Goff  and  Jordan,  1990],  an  algorithm  was  presented  for  stochastically  interpolating  2-D 
badiymettic  data.  In  this  report  we  specialize  to  1-D  bathymetric  profiles.  There  are  two 
motivations  for  this:  (1)  the  current  state  of  the  art  in  acoustic  modeling  is  best  developed 
for  1-D  and  (2)  we  can  generate  long,  highly  resolved  interpolations  without  the  extensive 
nesting  geometry  that  was  necessary,  due  to  space  limitations,  in  2-D. 

To  take  advantage  of  the  ARSRP  acoustic  teccmnaissance  experiment  we  have  selected  a 
wide  beam  profile  from  the  RV  Verna  cruise  2602  0/2602)  which  passes  through  the 
acoustic  reveiteraticMi  corridor  (Figure  1).  This  track  line  passes  very  near  three  of  the 
acoustic  target  sites  which  were  chosen  in  the  planning  stage  of  the  cruise  (Figure  1). 
Subsets  of  this  profile,  each  20  km  long,  representing  the  nearest  passage  of  the  V2602 
track  line  to  these  target  sites  (see  profile.  Figure  2)  were  chosen  for  stochastic 
interpolation. 

The  interpolated,  stochastic  component  of  the  profile  is  obtained  as  a  single  (arbitrary) 
realization  of  the  three-parameter  Gaussian-random-field  model  developed  by  Goff  and 
Jordan  [1988, 1989a]  to  represent  1-D  abyssal-hill  topography.  The  parameters  include  the 
rms  height,  H,  the  scale  parameter,  or  comer  wavenumber  k^,  and  the  fractal  dimension  D 
(valid  at  wavenumbers  much  greater  than  the  comer  wavenumber).  In  Section  3,  we 
review  the  stochastic  model  and  describe  some  of  its  mathematical  properties.  Numerical 
tests  of  the  inversion  algorithm  are  discussed  by  Goff  and  Jordan  [1989b],  and 
coDoparismis  between  synthetic  realizations  of  the  Gaussian  model  and  actual  Sea  Beam  and 
SeaMARC  n  data  sets  are  presented  by  Goff  and  Jordan  [1988]  and  Goff  et  al.  [1991]. 
The  data-syntitetic  ccnnparisons  indicate  that  the  model  is  generally  successful  at  matching 
the  low-order  characteristics  of  abyssal  hill  terrain.  Improvements  in  the  model  are 
anticipated,  primarily  through  improved  parameterizations  of  the  two-point  correlation 
frmction  and  the  addition  of  higher-order  terms  in  the  stochastic  characterization  [Goff, 
1990].  However,  it  is  hoped  that  the  preliminary  Gaussian  models  presented  here  will 
prove  useful  in  testing  numerical  schemes  for  modeling  acoustic  reverberations  generated 
by  small-scale  seafioor  topography. 
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Realizaticms  of  the  stochastic  model  are  used  to  interpolate,  to  1  m  resolution,  the  three 
20  km  subsets  chosen  from  the  V2602  track  line.  The  parameters  assumed  in  computing 
the  stochastic  realizations  ate  normally  estimated  by  applying  the  inversion  procedures  of 
and  Jordan  [1988]  to  multi-narrow  beam  swaths.  In  this  case,  however,  the 
resolution  of  the  single-wide  beam  data  is  insufficient  to  estimate  the  fractal  dimension,  and 
we  do  not  have  access  to  multibeam  data  in  the  vicinity.  For  the  realizations  generated 
here,  we  assume  a  fractal  dimension  based  on  other  estimates  in  the  North  Atlantic  [Goff, 
1991]  and  then  eyeball  the  rms  height  and  scale  parameter  for  which  the  covariance  model 
best  fits  the  covariance  estimated  from  the  V2062  profile.  In^rpolation  is  accomplished  by 
'molding'  the  synthetic  topography  to  the  constraints  provided  by  the  deterministic  profile, 
as  introduced  by  Goff  and  Jordan  [1990]  for  two  dimensional  fields.  The  net  effect  of 
molding  is  to  replace  the  smoother  features  of  the  unconstrained  stochastic  realization  with 
those  obtained  from  the  measured  profile.  The  molding  algorithm  is  described  in  Section 
4. 

One  of  the  sites  chosen  for  stochastic  interpolation  is  a  sediment  pond  (site  A  on  Figure 
2).  A  seismic  record  from  the  V2602  cruise  (Hgure  3)  reveals  the  large  scale  structure  of 
the  basement  of  this  pond.  In  this  case,  this  basement  structure  was  digitized  and 
interpolated  stochastically,  while  the  bathymetric  profile  (i.e.,  the  nearly  flat  sediment 
pond)  was  interpolated  linearly.  Interpolated  basement  bathymetry  and  sediment  cover  are 
fried  separately  so  that  their  geophysical  properties  can  be  distinguished.  No  substantial 
sediment  cover  can  be  detected  over  sites  B  and  C  in  the  seismic  profile  (Figure  3). 
However,  it  is  likely  that  a  thin  veneer  of  sediments,  undetectable  because  of  its  thinness, 
fills  in  the  nodes  and  cranies.  Filling  in  the  nooks  and  crannies  with  sediment  may  or  may 
not  prove  important  for  the  acoustic  models. 

To  model  the  possible  contribution  of  sediments  we  have  employed  an  efficient 
numerical  algorithm  for  modeling  the  sediment  distribution  on  very  rough  (fractal)  models 
of  seafloor  topography.  The  algorithm,  developed  by  Webb  and  Jordan  [1991, 1992],  is 
based  upem  a  cellular-automaton  description  of  sediment-particle  dynamics.  The  algorithm 
is  outlined  in  Section  S. 

Technical  specifications  regarding  the  V2602  profile  and  interpolated  subsets  are  given 
in  Section  6. 

2.  DETERMINISTIC  AND  STOCHASTIC  COMPONENTS  OF  SEAFLOOR  TOPOGRAPHY 

We  let  z(x)  be  the  height  of  the  seafloor  above  some  mean  reference  level  at  a  position  x, 
and  we  suppose  we  have  a  map  of  this  topography,  denoted  Because  the  map  is 
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based  on  limited  data,  it  may  accurately  represent  age-dependent  subsidence,  thermal 
swells,  major  fracture  zones,  oceanic  plateaus,  and  other  "large-scale"  features,  but  does 
not  ccmtain  topographic  variations  with  horizontal  dimensions  below  some  “cutoff  scale" 
Xfg.  We  assume  the  nu^  can  be  approximated  as  the  ouq)ut  of  some  filter  Af[z(x)].  To  the 
extent  that  the  mapping  cutoff  is  sharp  —  i.e..  M  passes  features  larger  than  with  no 
distortion  but  crmipletely  annihilates  features  smaller  than  —  this  filter  is  a  projection 
operator  z^i^x) »  M[zai(x)].  An  example  of  Zj^(x),  the  one  used  in  this  report,  is  the 
V26Q2  profile  which  has  been  gridded  at,  and  has  an  effective  cutoff  scale  of  1  km. 

We  seek  to  supplement  this  deterministic  description  of  the  seafloor  widi  some  stochastic 
representaticm  of  the  small-scale  features.  Let  h(x)  be  a  stochastic  process,  or  random 
field,  which  represents  the  statistics  of  the  topographic  variation  at  all  length  scales.  We 
define  A^fx) «  M[h(x)]  and  take  as  our  iiKxlei  of  seafloor  topography 

z(x)  =  ZA#(x)  +  h(x)  -  huix)  (1) 

In  other  words,  we  replace  the  stochastic  components  of  die  field  with  scale  lengths  greater 
than  the  cutoff  x^  by  the  known  ("deterministic")  components.  If  M  is  a  projection 
operatcff,  then  {^plying  it  to  this  model  recovers  the  map:  (x)]  =  z^ix). 

3.  ONE-DIMENSIONAL  GAUSSIAN  MODEL  OF  SMALL-SCALE  TOPOGRAPHY 

The  stochastic  model  used  in  this  report  is  a  stationary  (spatially  homogeneous) 
Gaussian  random  field  with  zero  mean  and  a  2-point  mcMnent  (covariance)  function, 

C«i(x)  =  {hixi)  h(xi  +  x))  (2) 

where  (-)  is  the  expected  value  and  x  is  the  spatial  lag  vector.  Under  the  Gaussian 
assun^tion,  all  higher  moments  of  the  random  field  can  be  expressed  in  terms  of  (2).  The 
power  spectrum  of  the  Gaussian  field  is  the  Fourier  transform  of  this  covariance  function 
[BracewelU  1978],  and  its  phase  spectrum  is  a  random  process  uniformly  distributed  on 
(0,2it]  [Priestly,  1981]. 

Although  the  distribution  of  seafloor  depths  often  fails  the  Kolmogorov-Smimov  test  for 
acceptance  of  the  Gaussian  hypothesis  [Gilbert  and  Malinverno,  1988],  approximating 
small-scale  topography  by  a  Gaussian  field  provides  a  simple  and  mathematically 
convment  description  of  its  most  important  features  [Bell,  197S;  Goff  and  Jordan,  1988]. 
These  include  its  RMS  height  variation,  the  orientation  and  characteristic  wavenumbers  of 
its  "tectonic  grain,"  and  the  variation  of  roughness  with  spatial  scale.  Moreover,  a 
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Gaussian  description  is  the  basis  for  the  study  of  higher-ofder  statistical  pr(q)erties,  which 
can  be  expressed  as  perturbations  from  a  Gaussian  form  [Goff,  1990, 1991]. 

Coworioitce  model.  The  mathematical  properties  the  two-dimensional  Gaussian  process 
enqdoyed  in  this  study  are  detailed  in  Goff  and  Jordan  [1988, 1989a;  1990].  This  model  is 
expressed  in  a  two  dimensional  form  to  account  for  azimuthal  anisotropy  (tectcmic  grain). 
Before  q)ecializing  to  1-D,  we  first  present  the  general  2-D  form  so  that  the  1-D  model  may 
be  placed  in  context  with  previous  work.  To  represent  smaU-scale  abyssal-hill  topography, 
we  have  proposed  a  covariance  function  of  the  frnm 

Caa(x)  =  /f2Cv(r(x))/Gv(0)  (3) 

where  H  is  the  rms  height  and  G^  is  defined  by 

Gyjr)  *  r^Kv(r),  0  S  r  <  «»,  v  e  [0,1]  (4) 

Ky  is  the  modified  Bessel  function  of  the  second  kind  of  order  v.  Gy  is  plotted  for  three 
values  of  vin  Hgure  4  (top  panel).  The  order  parameter  vcontrols  the  behavior  of  Gy(r)  at 
the  origin;  its  slope  at  r  *  0  is  zero  for  v  =  1  and  infinite  for  v  *  0.  Gi/2(r)  is  simply  an 
eiqxxiential  function. 

Azimuthal  variation  is  expressed  by  the  dimensionless  ellipsoidal  (Riemannian)  norm 

r(x)  =  [xT  Q  x]^'^  =  V<7nJC?  +  2qi2XiX2  +  qiixi  (5) 

The  scale  matrix,  Q,  can  be  expressed  in  terms  of  its  ordered  eigenvalues  ^  and  its 
normalized  eigenvectors  e„  and  , 

Q  e„e J  e,e  J  (6) 

Q  specifies  the  "outer  scale"  of  the  topography  through  the  characteristic  wavenumbers  k„ 
and  kg,  yielding  an  aspect  ratio  a  =  k^kg  for  the  lineation  or  "tectonic  grain"  of  the  abyssal 
hills.  The  orientation  of  this  anisotropy  is  given  by  the  directions  of  the  principal  axes. 
Since  the  covariance  decays  least  rapidly  along  the  e,  axis,  the  structure  tends  to  be  lineated 
in  this  direction. 

Thus,  five  parameters  determine  the  2-D  stochastic  model:  the  rms  height  H,  the  order 
parameter  v,  the  characteristic  wavenumbers  and  kg,  and  the  azimuth  Cs  of  in  the 
geographic  reference  frame.  The  scale  parameters  define  an  aspect  ratio  a  =  kjky 
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To  qwdalize  to  one  dimensxHial  surfaces,  we  choose  a  profile,  or  ship  track  directkm  €9 
where  dis  the  angle  between  the  and  ig directions.  The  lag  variable  r  is  then  exjmssed 
as  a  function  of  lag  jc  in  the  shq)  track  diiecdm: 

lix)  «  ksx,  kg  «  ^k}cos^0  +  ki  sin^a  «  [ejQeeP^  (7) 

Only  three  parameters  determine  the  1-D  stochastic  model:  rms  height  H,  the  cmler 
parameter  v,  and  the  scale  parameter  in  die  ship  track  directitm  k0. 

Power  spectrum.  The  1-D  power  spectrum  obtained  from  Fourier  transforming  (3)  with 
the  1-D  fonn  of  r  given  by  (7)  is 


PM)  -  ny  1-1/2) 

m  ,8, 

The  (Mie-dimensional  forms  of  the  power  spectrum  at  three  values  of  v  are  shown  in  Figure 
4  (bottom  panel).  Equation  (8)  is  a  power-law  spectrum  with  a  corner  wavenumber  At 
high  wavenumbers  the  power  spectrum  decays  at  a  rate  at  low  wavenumbers,  it 

is  flat 

Hausdorff  dimension  and  se^'Cffine  scaling.  The  Hausdorff  (fractal)  dimension  0  of  a 
topognqihic  surface  can  be  related  to  the  asymptotic  properties  of  the  covariance  function  at 
small  lag  [Adler,  1981].  Goff  and  Jordan  [19®^^  show  that  the  Hausdorff  dimension 
associated  with  (3)  is 


D  =  2  -  V  (9) 

Decreasing  the  parameter  v  increases  the  roughness,  with  the  limiting  cases  of  unity  and 
zero  corresponding  to  a  Euclidean  random  field  with  continuous  derivative  (D  =  1)  and  one 
winch  is  "space-filling"  (D  =  2),  respectively.  The  interum  case  D  =  l.S  (exponential 
covariance  function)  corresponds  to  a  Markov  process  (the  Qmstein-Uhlenbeck  process) 
[Fewer,  1971]. 

The  Hausdorff  dimension  describes  a  self-affine  scaling  relationship  at  wavenumbers 
much  larger  than  ke  [Goff  and  Jordan,  1988,  1989a].  A  topographic  surface  h{x)  is  seif- 
affine  if  there  exists  an  a  €  [0,1]  such  that,  for  all  F  >  0  the  topographic  difference 
function  d(x  -  Xq)  =  h{x)  -  h(xo)  is  identical  in  distribution  to  R‘*^diRx  -  Rxq).  For  the 


Goff,  webb,  and  Jordan 


ARSRP  Technical  report  2 


Page  7 


covariance  model  (3),  Goff  and  Jordan  [1988]  demonstrate  that  a=v.  The  self-afHne 
property  allows  the  topography  at  high  wavenumber  to  be  simply  interpolated  to  smaller 
scales.  However,  the  validity  of  this  interpolation  at  scales  below  the  resolution  of 
bathymetric  mapping  devices  is  hypothetical:  some  work  suggests  that  the  spectral 
exponent  [Fox  and  Hayes,  1985]  and  aspect  ratio  [Goffet  al.,  1991]  are  scale-variable.  On 
the  other  hand  analysis  of  a  deep-tow  record  near  the  Galapagos  ridge  indicates  that,  in  a 
rough,  recently  formed  terrain,  the  topography  exhibits  very  uniform  fractal  behavior  from 
the  kilometer  scale  down  to  the  10  meter  scale  [Goff  and  Kleinrock,  1991].  Small-scale 
stochastic  interpolation  using  the  covariance  tiKxlel  (3),  including  the  high-resolution 
synthetic  realizations  accompanying  this  technical  report,  require  testing  in  the  natural- 
laboramry  settings  that  will  be  investigated  during  the  ARSRP. 

Characteristic  scales.  Unlike  spectral  models  usually  associated  with  fractals  [e.g. 
Mandelbrot,  1983],  a  random  field  whose  second-order  properties  are  described  by  (3)  or 
(7)  does  not  have  infinite  power  at  zero  wavenumber.  Rather,  the  low-wavenumber  part  of 
the  spectrum  is  governed  by  a  characteristic  length,  or  outer  scale.  As  discussed  in  Goff 
and  Jordan  [1988],  the  characteristic  length  in  the  0  direction,  Xq,  can  be  defined  in  terms 
of  the  second  moment  of  the  covariance  function  in  the  9  direction,  which  yields 

,  2V  2(v  +  1/2) 

(10) 

where  is  the  scale  parameter  in  the  B  direction.  X^  is  interpreted  as  the  characteristic 
abyssal-hill  width,  and  Xj  as  the  characteristic  abyssal-hill  length. 

4.  MOLDING  ALGORITHM 

The  algorithm  for  stochastic  interpolation  involves  two  basic  steps:  (1)  generating  a 
synthetic  realization  at  the  desired  profile  length  and  data  spacing  using  the  stochastic 
model,  and  (2)  'molding'  the  synthetic  realization  to  the  constraints  provided  by  the 
deterministic  bathymetry.  To  generate  synthetic  realizations,  we  compute  the  Fourier 
spectrum  on  a  regularly  spaced  wavenumber  grid  by  multiplying  the  square  root  of  the 
power  spectrum  (8)  by  a  phase  factor  exp(i^),  where  ^  is  a  random  number  uniformly 
distributed  on  the  interval  [0,  27c)  [Priestly,  1981].  The  space  domain  image  is  then 
obtained  from  a  fast  Fourier  transform.  A  stochastic  realization  is  generated  such  that  its 
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total  length  and  sample  spacing  correspond  respectively  to  the  length  of  bathymetric  to  be 
interpolated  and  the  desired  resolution.  This  yields  the  /i(z)  term  in  equation  (1 ). 

The  nwlding  algorithm  contains  the  following  steps: 

1 .  The  bathymetiic  profile  is  regridded  at  a  regular  interval  (in  this  case,  1  km  spacing). 
The  portion  to  be  stochastically  intopolated  is  first  linearly  interpolated  (another 
inteiix>lating  scheme  may  be  used)  to  tiie  desired  resolution  (in  this  case,  1  m 
spacing).  This  step  yields  the  ZAf(jr)  term  of  equation  (1);  i.e.  it  represents  a 
raping  of  the  seafloor  depth  zCx). 

2.  The  stochastic  realization  is  subsampled  at  the  same  resolution  as  the  regridded 
bathymetric  profile  (1  km  spacing).  This  subsampled  realization  is  then  Unearly 
interpolated  (or  by  the  same  scheme  that  was  used  to  interpolate  the  bathymetric 
profile)  to  the  original  data  pacing  of  the  stochastic  realization  (1  m).  This  step 
yields  the  /tAf(x)  term  in  equation  (1). 

3.  hM(X)  is  subtracted  from  h(x)',  i.e.  the  large-scale  stochastic  morphology  is  stripped 
from  the  stochastic  realization.  What  remains  is  added  to  zm(x)  to  complete  the 
molding,  yielding  the  z(x)  of  equation  (1).  Note  that  z  (x)  =  zmM  where  x  is  equal 
to  one  of  the  original  (1  Ion)  grid  locations  (i.e.,  A/[z(x)]  =  zuix)). 


i 


i 
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The  following  model  parameters  were  used  to  generate  stochastic  realizations:  H  =  235 
nuke’=  0.45  km*l,  and  v  =  0.85  (D  =  1.15,  Ag  =  7.3  km).  Because  of  the  poor  resolution 
of  the  wide-beam  bathymetry,  v  could  not  directly  be  estimated  from  the  V2602  profile. 
However,  the  value  of  v  s  0.85  is  consistent  with  estimates  from  narrow  beam  records 
from  other  near-ridge,  unsedimented  sites  in  the  Atlantic  [Goff,  1990;  1991].  For  this 
value  of  V,  we  estimated  values  for  H  and  ke  simply  by  a  fit-by-eye  match  to  the  estimated 
covariance  for  a  portion  of  the  profile  in  the  vicinity  of  the  three  interpolated  sites  (Figure 
5). 

Figure  6  shows  the  results  of  the  molding  algorithm  applied  to  the  basement  morphology 
of  site  A  (compare  with  Figures  2  and  3).  In  this  figure,  the  dashed  line  represents  the 
topography  linearly  interpolated  from  the  1  km  gridded  wide-beam  bathymetry,  while  the 
solid  line  represents  the  1  m  stochastic  interpolation.  Note  that  at  the  location  of  the  1  km 
grid  points  both  profiles  coincide.  The  linearly  interpolated  sediment  cover  for  site  A, 
estimated  from  the  seismic  profrle  (Figure  3),  is  also  shown  in  Figure  6. 

5.  Sedimentation  ALGORITHM 

The  existence  of  sediments  may  have  an  important  effect  on  the  acoustic  response  of  the 
seafloor.  Sedimentation  rates  in  the  vicinity  of  the  acoustic  reverberation  corridor  are  on 
the  order  of  30  m/my  [Usitzin,  1972],  and  the  half-rate  of  seafloor  spreading  at  this  latitude 
is  ~15  km/my  [Tucholke  and  Schouten,  1988].  In  the  absence  of  any  reworking  of  the 
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sediment  cover,  we  would  expect  to  sec  -270  meters  sediment  cover  at  site  A  (-9  my 
crustal  age),  -510  meters  sediment  cover  at  site  B  (-17  my  crustal  age),  and  well  over  a 
kilometer  of  sediment  cover  at  the  oldest  crustal  ages  (-40  my)  spanned  by  the  acoustic 
reverbefation  corridor.  A  cursory  examinadon  of  the  V2602  seismic  profile  record  (Figure 
3)  reveals  a  much  more  ccxnplex  relationship  between  basement  age  and  the  sedimentation 
volume  in  this  region.  The  basin  containing  by  the  greatest  vdume  of  sediment  (up  to  500 
m  diick)  is  located  at  some  of  the  youngest  crustal  ages  on  this  record  (site  A).  The  paucity 
of  sediments  at  older  crustal  ages  (to  the  right  side  of  Figure  3)  and  the  fact  that  all  detected 
sediment  volumes  are  flat-lying  indicate  that  vigorous  post-depositional  sediment  transpon 
processes  are  operating.  No  sediments  can  be  detected  at  sites  B  and  C  on  the  1  km  scale 
(Figure  3),  but  we  can  assume  that  sediments  have  been  deposited  in  the  nooks  and 
oannies  at  the  1  m  scale  and  larger,  significantly  smoothing  the  surface  profile  at  those 
scales. 

Without  a  high-resolution,  two-dimensional  map  of  the  seafloor  topography  and  a  better 
of  understanding  of  the  nature  of  the  sediment  transport  processes,  it  is  difficult  to  quantify 
in  any  precise  way  the  distribution  of  pelagic  sediments  in  this  particular  area.  For  the 
purposes  of  prescribing  plausible  models  for  numerical  modeling  of  acoustic 
backscattering,  however,  it  is  possible  to  do  some  preliminary  calculations.  Webb  and 
Jordan  [1991, 1992]  have  developed  an  efficient  numerical  algorithm  for  calculating  the 
sediment  distribution  on  very  tough,  one-dimensional  models  of  seafloor  topography.  The 
algcuithm  is  based  upon  a  cellular-automation  description  of  sediment-particle  kinematics. 
It  consists  of  the  following  steps: 

1 .  Initialization.  A  specified  basement  topography  is  discretized  on  a  rectangular 
grid  with  vertical  spacing  &i,  horizontal  spacing  &,  and  total  horizontal  length  X  = 
KSx. 

2.  Sediment  Accumulation.  Parcels  of  dimensions  {Sx,Sh)  are  accumulated  at  each 

horizontal  position  xj^  for  M  time  steps,  yielding  a  topography  MSh. 

3.  Sediment  Redistribution.  Parcels  are  redistributed  for  2N  time  steps.  Transpon  during 
each  time  step  alternates  between  right-directed  and  left-directed  sediment  flow: 

a.  Right-directed  transport.  Beginning  at  k  =  1  and  moving  right  to  k  =  AT-l ,  is 

reset  to  where  =  min((hj^-h^^®^),  -min[(h^^j-/i^)/2,  0]). 

b.  L^t-directed  transport.  Beginning  atk^K  and  moving  left  to  /k  =  2,  is  reset 

U)  where  =  min((h^-hjj^®^),  -min[(/i^_j-hp/2,  0]}. 
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Tbe  sediment-water  interface  at  the  end  of  2N  time  steps  defines  a  topographic  profile 
which  remains  fixed  during  subsequent  sediment  accumulation  and 

T^stribution. 

4.  Steps  (2)  and  (3)  are  repeated  with  replacing 

5.  The  process  is  iterated  to  obtain  L  sediment  layers. 

Webb  and  Jordan  [1991, 1992]  have  demonstrated  that  this  cellular-automota  algorithm 
yields  a  diffusive  redistribution  process.  That  is,  sediments  are  transponed  across  a  flat 
hasonent  according  to  a  diffusion  equation  in  which  the  forcing  term  for  sediment  flux  is 
the  gradient  in  sediment  thickness.  This  transpon  is  characterized  by  a  diffusivity  ic  which 
scales  as  the  ratio  of  the  redistribution  time  to  ti»  accumulation  time, 

k-N/M  (11) 

For  a  small  value  of  k,  the  redistribution  rate  is  low,  and  sediments  drape  over  the 
basement  topography.  For  a  large  value,  the  rate  is  high,  and  sediments  are  swept 
downslt^  and  accumulate  in  topographic  lows  to  form  ponds. 

The  cellular-automota  algorithm  has  been  used  to  redistribute  sedin^nts  on  the  basement 
lopogitqihy  generated  by  stochastic  interpolation  far  sites  B  (Figure  7)  and  C  (Figure  8).  A 
seoiinent  volume  equal  to  a  layer  with  a  uniform  thickness  of  20  meters  was  assumed  for 
both  sites.  We  set  v  equal  to  4.7  x  10*3  m^/yr  for  site  B  and  2.0  x  10*3  rifi/yr  for  site  C. 
AldioQgh  no  diffusivities  have  been  estimated  for  pelagic  sedimentation  deep  oceanic 
settings,  these  values  are  on  the  order  of  those  observed  for  the  Bonneville  shoreline  scarp, 
a  region  of  semiarid  hillslope  degradation  [Cobnan  and  Waison,  1983;  Hanks  et.  al  1984]. 
They  are  sufficiently  large  that  sediment  ponding  occurs,  in  qualitative  agreement  with  the 
distribution  observed  on  high-resolution  proflles  [W.  Ryan,  personal  communication, 
1991]. 

It  should  be  noted  that  the  sediment  distributions  obtained  in  Figures  6  and  7  have  been 
generated  by  an  algorithm  that  allows  lateral  transport  in  only  one  dimension.  Allowing 
transport  down  topographic  slopes  in  the  cross-profile  direction  will  have  the  effect  of 
truddng  the  sediment  distribution  less  uniform.  It  is  possible  to  extend  the  Webb-Jordan 
algorithm  to  this  full  three-dimensional  problem  and  to  apply  it  to  basement  topography 
stochastically  extrapolated  out  of  the  plane  of  the  profile.  This  will  be  the  subject  of  a 
future  report. 
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6.  TECHNICAL  SPECmCATIONS  FOR  STOCHASTIC  INTERPOLATIONS 

The  following  model  parameters  were  used  to  generate  stochastic  realizations:  H  =  23S 
m, ikos 0.45  km'^  and  vs  0.85  (D  s  1.15,  7.3  km). 

The  fcdlowing  files  are  available  from  the  first  author,  all  in  ASCII  format: 

1 .  Original  V2602  profile,  'v2602.11d‘.  Header:  number  data  points.  Data  format: 
tetitiide  (decimal  degrees  N).  longitude  (decimal  degrees  E),  depth  (meters). 

(For  tite  remaining  files,  the  header  contains  the  following:  number  of  data  points,  data 
spacing  in  kilometers.  The  data  format  is:  depth  (meters).) 

2.  Regridded  V2602  profile,  'v2602.prof . 

3.  Stochastically  interpolated  site  A  basement, 'v2602a_si.base.  Site  A  corresponds  to 
points  indexed  1 1 1-131  of  regridded  V2602  profile. 

4.  Sediment  thickness  for  site  A,  'v2602a_si.sed'. 

5.  Stochastically  interpolated  site  B  basement,  'v2602b_si.base'.  Site  B  corresponds 
to  points  indexed  235-255  of  regridded  V2602  profile. 

6.  Sediment  thickness  for  site  B,  'v2602b_si.sed'. 

7.  Stochastically  interpolated  site  C  basement, ‘v2602c_si.base'.  Site  C  corresponds  to 
points  indexed  131-151  of  regridded  V2602  profile. 

8.  Sediment  thickness  for  site  C,  'v2602c_si.sed'. 

The  surface  profile  for  any  site  is  generated  simply  by  adding  the  '.base'  and  '.sed' 
files.  To  obtain  files  either  by  telnet  mail  or  FTP,  contact  John  Goff  at 
goffi@jeffreys.whoi.edu  (telnet  address  128.128.16.216). 
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8.  FIGURES 


Figure  1.  Track  line  of  the  V2602  bathynKtric  profile  and  seismic  record  through  the 
ARSRP  acoustic  reverberation  corridor  (dashed  box).  Tick  marks  are  spaced  every  1(X) 
data  points,  x's  within  box  represent  preliminary  acoustic  target  sites  for  the  upcomming 
ARSRP  acoustic  reconnaissance  experiment.  The  Mid- Atlantic  Ridge  plate  boudary  is  also 
shown,  with  double  lines  representing  spreading  ridge  segments,  and  single  lines 
representing  ridge  offsets. 

Figure  2.  V2602  wide  beam  bathymetric  profile.  For  consistency  with  the  seismic 
proHle  (Figure  3),  East  is  to  the  left.  20  km  sections  chosen  for  stochastic  interpolation  are 
marked  as  sites  A,  B,  and  C.  These  sites  correspond  to  the  nearest  points  of  passage  to  the 
respective  sites  A,  B,  and  C  marked  in  Figure  1.  The  axis  of  the  Mid-Atlantic  Ridge 
(MAR)  is  also  marked. 

Figure  3.  V2602  seismic  profile.  Compare  with  bathymetric  profile  in  Figure  2.  This 
record  was  used  to  digitize  the  basement  depths  for  site  A.  Note  general  absence  of 
sediments  at  older  crustal  ages  (crustal  age  progresses  left-to-right). 

Figure  4.  Functional  form  of  the  model  covariance  function  Cv(r),  plotted  for  values  of 
V  =  0, 1/2,  and  1  (top  panel),  and  their  normalized  Fourier  transforms  plotted  on  a  log-log 
scale  (bottom  panel). 

Rgure  5.  Covariance  function  estimated  from  a  200  km  section  of  the  V2602  profile 
(solid  curve)  and  fit-by-eye  covariance  model  (dashed  curve)  (Equation  3).  The  200  km 
subset  begins  at  140  km  and  ends  at  340  km  from  the  beginning  of  the  V2602  profile.  The 
parameters  used  in  the  covariance  model  are  //  =  195  m,  ke  =  0.45  km'l,  and  v=  0.85  (D 
=  1.15,  Xq  =  7.3  km). 


Figure  6.  Stochastically  interpolated  basement  depth  for  site  A  (bottom  solid  curve)  with 
measured  sediment  cover  (top  solid  curve).  Dashed  lines  represent  1-km  gridded  basement 
dq)ths  which  were  used  to  constrain  the  large-scale  behavior  of  the  stochastic  realization. 

Figure  7.  Stochastically  interpolated  basement  depth  for  site  B  (bottom  solid  cuvre)  with 
numerically  generated  sediment  cover  (top  solid  curve).  The  average  sediment  cover  is  20 
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m.  Dashed  lines  represent  1-km  giidded  basement  depths  which  were  used  to  constrain  the 
large-scale  behavior  of  the  stochastic  realization. 

Figure  8.  Stochastically  interpolated  basement  depth  for  site  C  (bottom  solid  curve)  with 
numerically  generated  sediment  cover  (top  solid  curve).  The  average  sediment  cover  is  20 
m.  Dashed  lines  represent  1-km  gridded  basement  depths  which  were  used  to  constrain  the 
large-scale  behaviOT  of  the  stochastic  realizaticx). 
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